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I. INTRODUCTION 


A. BACKGROUND 

Elastic wave propogation in a fluid-saturated porous solid 1s a very 
important and active subject of research with a wide spectrum of applications 
including oil prospecting, nuclear test ban treaty verification and sound 
transmission and reflection at the ocean floor. 

When the fluid in the pores’ is highly compressible (air for example) the 
acoustic properties are completely determined by the characteristic impedance, 
a - or, equivalently by the complex wave number, k . From these quantities 
can be extracted parameters which depend only on the microscopic geometry of 
the porous solid, the values of which can then be used when the fluid in the 
pores 1s not highly compressible (water for example). It is therefore 
desirable to have accurate and reliable methods by which these parameters can 
be measured. 

Be HISTORY OF CHARACTERISTIC IMPEDANCE AND COMPLEX WAVE NUMBER MEASURING 
mECHNIOQUES 

Rayleigh (1929) [Ref. 1, p. 328] has investigated the performance of an 
idealized absorber of sound consisting of small bundles of tubes. He extended 
the theories of Helmholtz and Kirchoff that relate to the transmission of 
sound along tubes and was able to deduce a value for the absorption 
coefficient at the surface of such a material in terms of the propagation of 
the surface area occupied buy the open pores, the radius of the tubes, and the 
physical constants of air. This analysis has provided a basis on which to 
found more general theories and description of sound propagation through 


porous materials. 


tk 


Scott (1945) [Ref. 2] pronosed a method of determining the complex wave 
number, & which Gescribes the spatial variation of the acoustic disturbance in 
a porous medium. This method consists of measuring the attenuation constant 
and the wavelength constant of a plane progressive wave by nrobing the porous 
medium with a probe tube microphone. Scott then determined the chatacteristic 
impedance of the material with an annaratus devised by H. O. Taylor (1913) 
[Ref. 3] to measure the sound damping effects of acoustical materials. This 


apparatus 1s similiar to the present day impedance tube. 





Beranek (1949) [Ref. 4] describes several methods of obtaining a | 


material's acoustic imnpedance by measuring the properties of the sound field 
in the vicinity of an absorbing material exposed to a normally incident sound 
wave. Two of the most tractable methods are based on transmission line 
theory: the impedance tube method and the resonance method (normal 
incidence), with the latter method generally being viable only for determining 
Surface impedance. Beranek also developed a method to calculate the complex 
propagation constant from measurements of a material's characteristic 
properties, 1.@€. porosity, pore size, etc. Once these properties are 
determined, the characteristic lmpedance can be calculated. 

More recently, Yaniv (1973) (Ref. 5] developed a new method of 
calculating both the complex propagation constant y = jk and characteristic 
impedance of a material using transmission line theory applied to an impedance 
tube device. His method consists of measuring standing wave properties in the 
proximity of an absorbing material alternately backed with a high impedance 
and a low impedance termination. 

Finally, in 1980, Chunq and Blaser [Ref. 6] published their article 
discussing a Dual-Microphone Transfer Function technique for measuring in-duct 


acoustic properties, namely the acoustic impendance and_ reflection 


12 


coefficient from which the complex wave number can be determined. This 
technique waS compared to the conventional Standing Wave Tube technique and 


found to give similar results [Ref. 7]. 


fa PURPOSE OF STUDY 

The objective of this research 1S to investigate the feasibility of 
obtaining the characteristic impedance and the complex wave number of a fluid- 
filled porous solid by measuring the resonance frequencies and half-power 
widths of a fluid-filled rectangular cavity which is partially occupied by a 
porous solid of varying thickness. In this investigation, the fluid is air, 
the porous solid is fiberglass and the cavity is made of plexiglass. 

This investigation differs from previous investigations in that 1) Biot 
Theory [Ref. 8, 9, 10, and 11] is applied to solve for the allowed modes of 
the partially-filled cavity, and 2) for the primary mode of interest the 
particle velocity 1S directed along the length of the fluid-porous solid 
boundary aS opposed to normal to this boundary. 

Table 1 sSunmarizes the advantages in the course of this research. The 
advantage of this technique over the currently existing methods is that it 
does not require probing the sound field either within or in front of the 
porous solid as do methods 1, 2, 3, and 5, nor does it require a pair of 
phase-matched microphones aS does method 6. Additionally, characteristic 
properties of an absorbing material need not be known as in method 4. Also, 
the acoustic properties of the bulk fiberglass can be determined using only 
thin samples. 

The theory developed for the new method is presented in Chapter II along 
with a general overview of the equations used in the classical methods. The 
experimental results for each of the methods used in comparison are presented 
in Chapter III and comparisons are made and discussed in Chapter IV. 


13 


TABLE 1 


ADVANTAGES AND DISADVANTAGES OF METHODS USED TO DETERMINE THE CHARACTERISTIC 


‘Method 


1. Probe Tube 
(Scott) 


2. Impedance Tube 
(Beranek ) 


3. Resonance 
(Beranek ) 


4. Material Properties 
(Beranek ) 


IMPEDANCE AND COMPLEX WAVE 

Advantages 

a) Simple theory. 

b) Simple procedure. 

a) Simple theory. 

b) Simple procedure. 

c) Probing of material 
not required. 

d) Normal absorption easy 
to measure. 

a) Probing not required. 

b) Measurement of re- 
sonance frequency and 
half power widths is 
precise. 

a) Sound wave not needed. 
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NUMBER 


a) 


a) 


b) 


c) 


a) 


b) 


Disadvantages 


Thick sample needed 
to get character- 
istic impedance 
without a wave 
mode 1 for the 
porous medium. 


Must 
field. 


probe sound 


Thick sample needed 
to get character- 


1SELE impedance 
without a wave 
model for porous 
medium. 

Theory 1S somewhat 
Gitt icu lee 

Discrete  frequen- 
cles. 

Measures surface 
impedance. 


Requires low losses 
for well defined 
resonances. 


Complex theory. 
Measurement of 


microscopic proper- 
ties 1s complex. 


Method 


Bee Yaniv 


6. Two Point Phase Dif- 
erence (Chunq and 
Blaser) 


7. New Method 
(Raker) 


AC 


(cont inued ) 


Advantages 


Simple theory. 


Simple procedure. 


Probing of material 


MOtwteeeoocadty .« 


Simple theory. 


Simple procedure. 


Probing of material 


not necessary. 


Probing not required. 


Simple procedure. 


Ability to determine 
DC flow resistance 
non-invasively. 


Bulk properties 
obtained from Ean 
samples. 


libs) 


Cc) 


Db) 


Disadvantages 


Requires classifi- 
Gatton Ob material 
AS mile Geom sot t.. 


Time consuming. 


Quarter wavelength 
termination needs 
continuous adjust- 
ment with frequen- 


Ey. 


Thick sample 
needed. 

Microphone posi1- 
tions only optinum 


for select frequen- 
cies. 


Requires accurate 
microphone phase 
calibrations. 
Complex theory and 
numerical methods 
required. 

Discrete frequen- 
eles. 


Requires low losses 
for well defined 


EeESOnaNCES. 


Tia Pa EOERs 


The theory behind the new method 1s presented in detail in this section 
with much of the fundamental concepts and understanding having come from the 
investigation of sound propagation in a porous. solid saturated with 
suner-fluid Helium [Ref 11]. Additionally, the basic theory and equations used 
are presented for each of the classical methods used in comparison with the 


new techniques. 


A. THEORY OF A ELUID-SATURATED RIGID, PoR@US 2orie 


The linearized equation of motion tor the microscopic fluid velocity 


ae ee in the Eulerian description is [Ref. 12, p. 49], 
* (2.15 
Vmicro a aie ee av 1 7( 7eV. 
Pt: eae + “nicro’Y)Ymicro} = 7% + MV nicro + (4 ay VTS nicro)s 
Using 
ne 2 > a > > S 
Vmicro = “V'’micro) 7 YX YX Mnicro | (2.2) 
we can write Eq. (2.1) as 
37 . > > > > > > > > > > 
micro ae ec, =o vo a “Va 
PEI 5E + (Ymicro* 9)" micr w+ nV x VX Yniorg + (6 + > nV VS 
(2338 


Where pe is the fluid mass density, p is the fluid pressure; 7 1Seei fluid 
shear viscosity and ¢ is the fluid bulk viscocity. The last term in this 
equation can be ingnored here because its affect 1s small compared to the 


other terms on the scale of a pore. Ignoring the last term, Eq. (2.3) becomes 


> 
DV > >> > 
MLSE ee = nVX VXV 


Zeal ag MiCce um (2.4) 
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> 
DVmicro . ; oe . 
where p> is the total time deriviative of the microscopic Beoiicl 


velocity. From Eq. (2.4), a phenomenological, linearized equation for the 
average (over the microscopic structure) fluid velocity, denoted < ey can 


be written: 


Gitioas Ran % 
iC liens cnn sn a ~V<p> = RelowS “micro” ; (2.5) 
Here perf is the phenomenological effective fluid mass density and R¢),, is 
the phenomenological flow resistance [Ref. 13, pp. 252-255]. Re¢jow iS defined 


such that in steady state 
> 
-¥«p> = Roc “micro? - (2.6) 


In general (including oscillatory flow) pe¢¢ and R—ej,4, can be operationally 


defined by 
> 2 
> Pett een = Kinetic ic opneray per unit (27) 
Kee ss = Rate of dissipation of (2.8) 
Rflow S“micro kinetic energy density 


For eJ4t time dependance, Eq. (2.5) can be written, upon combining both 


velocity terms on the left hand side, 


= 
E pepe) + Rerow! | “micro? =~ YP? + (2.9) 
where pe¢f and RejQ, are now frequency dependent phenomenological properties 
which describe the effect on the fluid due to the microscopic geometry ([Ref. 
Il, pp. 334-335]. It should be noted that in Eq. (2.9) and in the remainder of 


this thesis complex quantities are denoted by a tilde. 


Ly 


Defining complex effective density, peg¢¢e(@), as 


~ ~ ~ R ( w) 
pef E(w) = Pofe(w) + — ) (2.10% 
Eq. (2.9) may be written in a more comnact form as 
aoe = = 2 
tears! w) SY ie cee = = /op> Se (2.18 


Plane progressive wave solutions to Eq. (2.11) travelling in the positive 
x-direction may be sought by letting the acoustic pressure be given by 
sp el(ut - kx) ang letting <Vniieror® vane ej(ut - Kx) | where k = k - Ja with 
k, a > 0 Such that the energy of the wave exponentially decays as _ one 
moves deener into the porous’ solid (note in order to avoid confusion later, 
the letter a will be used to denote the attenuation constant, reserving 


the letter a for the tortuosity in the Biot theory). Eq. (2.11) then 


becomes 


JO PeglO) Vi =i Kap . (2.129 

The microscopic equation of continuity for “micro Still holds for vs oral 
if we consider a large enough volume so that irreqularities on the order of 
the pore size average out, and if we consider the density of the fluid in the 
pores to oscillate about the bulk fluid density with the passing of an 


acoustic disturbance. Hence we let 6 = p¢ + sp elf ut - Kx) and the equation 


of continuity becomes, to first order, 


nt 


jwsp - jpg Kk Vayqg = O (2.138 


g 


< 


The relation between pressure and density variations (the equation of 
State) in a fluid-saturated porous solid can be written 


2 (uw) 


Ey 





} (2.148 
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where C¢r 1s the speed of sound in the bulk fluid, ay is the complex ratio 
of specific heats (cp/cy,) as a function of frequency, and y is the adiabatic 
ratio of specific heats. 

Eq's. (2.12), (2.13) and (2.14) lead to a wave equation with the phase 


velocity determined as follows: writing Eq. (2.12) (the equation of motion) 





as 
ef eo) Vavg = 2 69 a) 
and substituting Eq. (2.14) into it results in 
: Pete ®) Yayg = Cg? ie) 6p OIG 


and dividing Eq. (2.16) 


LOM. 


Recasting Eq. (2.13) (continuity) as p¢ Vavg oe § 
by it results in 





w Oo 
(pe)? = ce? = - Keel (2.17) 
Pf 


Defining the complex tortuosity in the fluid saturated porous solid by [Ref. 
pee. 135 and Ref. il, p. 37] 


~~ Bagel) 
eo (2.18) 


Pf 


the phase velocity can be written as 





yo cf Wa) 
oa? = W)\? = G) 
The complex tortuosity can be written as 
ere Rak ( o) 
oleae (2.20) 
J WPF 
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where a is the intrinsic, or infinite frequency (inviscid) tortuosity, and 


a id 


F(w) is a dimensionless correction function whose purpose is to make the 
frequency dependence of a(w) come out right [Ref. 8, p. 185 and Ref. 11, p. 
BT lw 


The complex tortuosity can be written in the form [Ref. ll, p. 49] 


a(w) = a- jrFlw) , (2.21 
where 
_ R 
=o (2.22) 
WPF 


is a Gimensionless flow resistance: 
The mathematical form for the dimensionless correction function Gane 


taken as [Ref. ll, p. 57 and Ref. 14) 


F(w) = 1+ iG)? (2.23) 


where 
: (2.253 


Here 6 1s the structural factor in Biot's universal correction function Vandi 
1s the intrinsic tortuosity of the microscopic pore geometry [Ref. 10]. 

The intrinsic tortuosity is a weighted path length correction factor which 
Satisfies a2l . It depends most strongly on porosity P , with al as 
P+l (as is the case for fiberglass). Typical values range from 1 to 4 [Ref. 
lye 227 sand Ret mls. 


The DC flow resistance is proportional to n/(pore size)*, so that 





ia Saal 1 Z (Viscous nenetration depth) 2 (2.25) 
a aopp ° (pore size) ¢ pore size ° i 
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The quantity r/o determines the relative importance of flow resistance versus 
mira inertia in determining the microscopic flow field. 

The structural § depends upon the average cross-sectional shane of a 
pore [Ref. 8, p. 185]. Experimentally it appears that 624 for many granular 
merous solids [Ref. 16]. 

eMecereaeimg &q. (2.22), (2.23) and (2.24) into Eq. (2.21) results in 


BNE eee 


Roe os oe ee oer 








a(w) = a{l-j ane 


moerall calculations in this research it will be assumed that a=l and 6=4 for 


fiberglass [Ref. ll, pmo. 49-59]. 


The complex ratio of specific heats, y(w) , is defined as 


¥(w) = Cp/Cy (2.27) 
where cy, and c, are, respectively, the complex, frequency dependent specific 
heats per gram at constant pressure and at constant volume. For an ideal gas 


~ ~ _R . 
ey where M is the molar mass, and so 


y(w) = 1+ (2.28) 


Mey 
The pronagation of sound in the porous medium will vary from isothermal 
at low frequencies when the thermal penetration depth is large comnared to 
a pore size, hence Cy +o, to adiabatic at high frequencies when the thermal 
penetration depth is small compared to a pore size, and Cy = bulk fluid 
value. Since the transport of momentum and energy in a gas both arise from 
molecular collisions, it will be assumed that thefunctional form 1s) - 1) has 
the same form as the complex tortuosity Pe only with the thermal 
penetration depth &}, in place of the viscous penetration depth 6,;.,. 


Hence we take [S. Baker, private communication] 


aan eae, 
ye = ol 1-3 FG) ee (2.29) 


ie 


Vigan, 


a 


rt 


where r and F(w) are given by Eq's (2.22) and” (2.23), Yrespectiveliyaumaaae 


viscous penetration depth, Syigc » and the thermal penetration depth §,, are 


2 2a z aes 
Ovisc = wer Sth = WPECP (2.308 


where n is the fluid shear viscosity, pg is the bulk fluid density, Kk “ig 


defined by 





the thermal conductivity and Ch is the specific heat per gram at constant 





oressure. 
Sth 
The factor —— Can be recast as 
Svisc 
Sth 7 eel 
; = SG Soe (2.3 
Svisc 
where 
v = r/pg = kinematic viscosity (2.329 
and 
x = K/p Ga thermal diffusivity (25338 


and Pr = Prandtl number which has a value of 0.733 for air at 20°C. Hence the 


complex ratio of specific heats can be rewritten as 


ieee (2.34) 





In the high frequency limit, the square of the phase velocity in the 
porous material (Eq. 2.19) qoes to the fluid inertia dominated limit in which 
sound propagates at the bulk fluid speed reduced by the square root of the 
tortuosity. Conversely, in the low frequency limit, flow resistance 


a) 





dominates, and pressure oscillations diffuse through the porous material. The 
crossover from inertia dominance to flow resistance dominance occurs when the 
viscous penetration depth is approximately equal to a characteristic pore 


size. Figure 2.1 summarizes the limiting conditions. 


1. Development of Dispersion Relation in a Fluid-Filled Cavity Partiall 

Filled With a Porous Solid 

Consider a rectangular cavity which is partially filled with a fluid- 
Saturated porous solid. The region occupied by bulk fluid is denoted as 
region 1; the region occupied by the fluid-saturated porous solid is 
identified as region 2. The porous solid is assumed homogeneous’ and 
PSOLropic. Figure 2.2 depicts the geometry used to develop the dispersion 
relation for standing waves in the cavity. It will be assumed that all cavity 
cavity walls are perfectly rigid so that neu = 0. Then Scie = 0 at these 
boundaries and the velocity potential in regions 1 and 2 can te written as 
follows: 


In region 1 (bulk fluid), 


J wt 


1 = $190 coslk x] cos[kyz(z-d1)] (2235) 
where 
eee = 
Ky rf n Eee 6. Soret (2.36) 
w2 = Cea kiz2) cE (oe7) 
and 


bulk modulus of fluid 


Be 


d;9 = complex amplitude of the velocity potential in region l. 


> 


High and Low Frequency Limits on Phase Velocity 
in a Porous Material 


+ Fluid Inertia Dominates 


High Frequency Limit Low Frequency Limit 
a = 
——— ie = Ce 
ap os kr ee * es 
= a eed —————— 
Ee eee 
mre Le 
See a 
ee eae a at 
ae ee eee 
a _ eae 
—S— 
Plug Flow Laminar Flow 
EB) > cee a( a) > 
Pf 
y(w) + y y(w) > 1 
2 aoe aoe ccee 


Flow Resistance Dominates 


+ Sound Propagates + Sound does not Propagate 
(Pure Diffusion) 


Cross~over when: 


Svyisc 7 Characteristic pore size 


Figure 2.1 High and Low Frequency Limits on Phase Velocity in a Porous 
Material 
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REGION 1 
(Q@ULK FLUID) 


tw > 


REGION 2 CPoROUS MATERIAL) 


Figure 2.2 Rectangular Cavity Geometry for New Method. 


ZS 


The x-component of the wave number 1s a pure real number since the houndary 
conditions of rigid walls normal to the x-direction at x = 0 and == 


real and discrete values for ky - 





In region 1 the narticle velocity can be written as 
> _ SAS ; 
vies ee (2.398 
and the pressure as 
S51 = ~jW O; $1 - (2.40) 


In region 2 (porous material), 








o2 = $20 € 3“ cos [kx] cos [k9z(z+d2)] (2.41) 
where 
apt -3 = 
kK. aad L f rn QO, iy Z5 ° 6 ° f (2.425 
oe Zz 22 
we = (ko + Koz Neg , (2.439 
2 
~ 2 CE y( w) 
and 


$20 = complex amplitude of the velocity potential in region 2. 


The particle velocity in region 2 is 


<S¥ 2 


2 = Von (2.45) 


ana the ‘pressure is, from Ba. (2-11), 


S52 = -jwalw) o¢ o9 . (2.46) 
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With the velocity potentials defined in both regions 1 and 2 all that 
is left is to match conditions at the boundary between the bulk fluid and 
porous material (i.e., at z = 0) to produce the desired dispersion relation. 
The two boundary conditions that must be satisfied at 2=0 are: 


1. continuity of pressure across the boundary: 


$,(z=o) = a uw) (z=0) (22473 


and 


2. continuity of the normal component of volume velocity: 


361 30 
= = p i : (2.48) 


=0 2-0 


where P is the porosity of the material. Applying boundary condition 1 


results in 

$19 coslkyzd,] = a(w) $29 cos{kzd2) . (2.49) 
Applying boundary condition 2 results in 

dio Kiz sin (kizdil = -$20 P Kaz sin[kgzd2] . (2.50) 


Dividing Eq. (2.50) by Eq. (2.49) results in the dispersion relation 





E 
Kiz tan [ky2d 1] ar pe coe tan ed) =O , (26515) 
where 
S fe ml 172 





my 
Koz = a a( w) a a kK ’ (2.59 
f 


and 
k es i= O; i 2 ee e¢ e e (2.545 


“te ON “te Ne 


In Eq. (2.53) the forms of a(w) and y(w) have already been defined in Eq. 
(2.26) and Eq. (2.34), respectively. Eq. (2.51) is a transcendental equation 
in é whose roots determine the natural frequencies and half power widths or, 
equivalently, the quality factors, of the allowed "standing wave" modes of 
this system. From these roots can be determined the complex wave numbers 
ae and oo and the resulting pressure distribution. 

Two approaches to the use of the dispersion relation become evident. 
The first 1s the "forward" problem, where the roots of Eq. (2.51) for a may 
be regarded as functions of Rpc and a , and measured resonance frequencies 
and O's may be compared with predictions; the second in the “inverse” 
problem, where Rpc and a are regarded as functions of the roots of Eq. 
(Zoi) Gror : and are found by variation such that the best fit is obtained 
between the measured and calculated spectra. The aim of this work 1s to 
verify the technique to solve the inverse problem by comparing the value of 
Roc extracted from resonance data with that actually measured. Appendix 
A discusses the numerical technique to determine the roots of the dispersion 
relation and the algorithm used to determine the best fit value of Ro: 

Once the "best fit" value of Roc has been determined, then, since 


the roots of the dispersion relation in terms of w are known for this value 


of Roc , the complex wave number is found from 


QUE 2 


k = 


Pie  EOUnc Ero ha'S..9 (2.19), (2.26) and (2.34). The characteristic 


impedance can be calculated from 


Zenee a Pore (”) é (2257) 


where Pace) is civen by Eq. (2.10). Noramlizing the characteristic 


impedance to that of the bulk fluid (foc) results in 


Zchar (normalized) = Zchar/POS ° (2.58) 
The complex pressure field for both regions may be calculated once 
the roots of Eq. (2.51) are found. The pressure distribution in region 1 can 


be written in terms of the velocity potential in reqion 1 as 
5, = ire A eee) 


where $1 is given by Eq. (2.35). The pressure distribution in region 2 can 


be written in terms of the velocity potential in region 2 as 


- a - 40 X(B) pe d, (2.60) 


where $5 is given by Eq. (2.41). 
The velocity potentials Dy and b> are related by the boundary 


condition that pressure be continuous’ across the interface between the bulk 


fluid and the porous material, Eq. (2.47), which requires that 
Dio cos [k, 24, ] = & (a) Doo cos[Ka7de]_ . (2-617) 


For simplicity, choose D9 to depend on $j) , resulting in 


d a Dio cos ky 2d, 


20“ €(B) cos Kaz d, aa 


The pressure distribution in each region may be written as follows: 


SS, = Ec Pe cos [kx] cos [k, Aah " Bio eJ°t (2263) 


Ze 


and for region 2, 


Cos kK d 2a 


a cos [k55(z+d,)] oe eos (2.64) | 


2a SID ee CoS eae 0249 


With Eq. (2.63) and Eq. (2.64), and the values of kj, amd ko, determined 
from the roots of the dispersion relation of Eq. (2.51), the pressure 
distribution throughout the cavity may be calculated. See Appendix 8B for a 


more detailed discussion of the pressure calculation. 


B. CURRENT METHODS USING TRANSMISSION LINE THE@GEY 

Three of the five classical methods for determining the characteristic 
impedance and complex wave number of a porous acoustic material are based on 
transmission line theory. These methods involve measurement of various 
properties of the standing wave and data are taken at various points somewhat 
removed from the sample surface [Ref. 4, pp. 317-321]. The three methods 
considered are: 

-- Impedance Tube Method 

-- Resonance Analysis Method 

-~ Yaniv's Method 
Only Yaniv's method allows direct determination of the material's 
characteristic impedance and complex wave number. The impedance tube method 
allows direct measurement of the resistive and reactive components of the 
surface impedance. The complex wave number can only be calculated if 
reflections from the far side of the porous material can be neglected and one 
has some knowledge of the complex effective mass density of the fluid in the 
pores. As developed by Beranek, the resonance anlaysis method measures the 


same quantities as the impedance tube method, with the same limitations. It is 
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discussed primarily for its similarity to the new technique in that the 
parameters of the resonance curve are measured. 
1. Impedance Tube Method 

In this method the non-dissipative acoustic transmission line is 
considered a smooth rigid-walled tube whose lateral dimensions are small 
compared to a wavelength. The losses at the wall and in the bulk volume are 
considered small compared to the losses associated with the absorbing material 
under test. Assuming that the source of sound is at one end of the tube (x=0) 
and the sample under test is located at the opposite end (x=L), the surface 
impedance can be determined. Figure 2.3 depicts the geometry of the impedance 
tube. 

The impedance tube that has a constant cross section S and_ the 
source vibrates harmonically at a frequency sufficiently low that only plane 
waves can propagate. A standing wave pattern is set up by the right and left 
travelling waves and can be expressed mathematically as 


~ _ @ itu + k(L-x)) , = Jifut - k(L-x)] 


p (205) 


od 


where 0 is the acoustic pressure, A and B are determined from the boundary 
conditions at x=0 and x=L. 

It has been shown [Ref. 17, pp. 200-206] that the surface impedance 
of the termination may be written in the form, 


Zt, Z 1+ B/A el 


a (2.66) 
Om) 8/5 ke 


where B/A is the ratio of the relected wave amplitude to the incident wave 
amplitude and 6 is the phase difference between the relected and incident 
wave. Thus given B/A and 6 , the surface impedance of the termination can be 
calculated. 


Bi 


Solving for the amplitude ID | of the standing wave it can be seen that 
the amplitude at pressure antinode is A+B, and at a pressude node is A-B. The 
ratio of the pressure amplitude at an antinode to that at a node is the 
Standing wave ratio, or SWR: 


+ Panti 
ALB ee Samer (2.67) 


UR 3-B ~ Dnode ’ 


which can be rearranged to give the ratio B/A: 


B _ SWR-1 
Neat SWR+] e (2.68) 





Thus measurement of the SWR by probing the sound field in the impedance tube 


with a small microphone yields the value of B/A [Ref. 17, p. 206]. 





Figure 2.3 Geometry for the Impedance Tube 


The phase angle 8 can be evaluated from determining the distance to 
the first node from the end at x=L. These nodes are located at 
k(L=x,) = 9/2 = (n-1/2)i Gon = ee (2.628 
so that the phase angle for the first node is 


@) = 2k (Lex decom (2.70) 


S2 


With the SWR and 6, known, the resistive and reactive portions of the 
impedance can be determined graphically using a "Smith" chart. To determine 
the characteristic impedance of the porous material, the sample thickness is 
increased such that the reflected wave from the rigid termination is 
insignificant compared to the incident wave. From the definition of the 
material's characteristic impedance, the complex wave number can _ be found 
assuming that we know how the fluid reacts in the pores of the material. The 


characteristic impedance may be written as 


~w ~~? ~w 


Zchar = Peff © (2.71) 


where baer is the complex effective fluid density in the pores and c is the 
phase velocity of the acoustic distrurbance in the pores. The phase speed 


C in a fluid, can expressed as 

c2 = Be/dece (2.72) 
where By is the complex bulk modulus of the fluid in the pores. Therefore 
Eq. (2.71) becomes, upon substituting Eq. (2.72) for Peer 


Zchar = Be/c ee) 


~w 


and since @ Se the equation that describes the complex wave number R can 


be written as 


k = B- Foha : (2.74) 
If is is assumed that the material pore size is small compared to the viscous 
and thermal penetration depths (in air these two parameters are approximately 
equal) in the frequency range of interest, we can write 

Be = By (2.75) 
where Br is the isothermal bulk modulus of the fluid in the pores. 


Additionally for an ideal gas, the isothermal bulk modulus can be written as 


S¢ 


Br = l/«p (2.76) 
where xp 1s the isothermal compressibility of the gas, and is equal to 1/P, 
where P is the pressure in the pores [Ref. 18, p. 47]. Hence Eq. (2.74) can 


be written, assuming isothermal conditions exist, in the form, 


k = SUSIE Zee 7 
where k =k - ja, for a wave travelling in the positive x-direction in the 
porous solid. Eq. (2.77) is presented to show the relationship between the 
complex wave number and the characteristic impedance assuming isothermal 
conditions in the pore volume. This equation was used only to check Beranek's 
Material Properties Method in the isothermal limit. Only the values for 


7 eee extracted by the impedance tube method are reported here. 


2. Resonance Analysis Method 


The resonance analysis method is used to determine the surface 
impedance of a porous acoustic material. This method requires low loss within 
the material to be successful. A resonance curve may be traced out by varying 
either the tube length or the frequency in this’ method. For the present 
research, the impedance tube discussed previously was used, with frequency 
swept to trace out a resonance. 

Beranek presents the equation needed to determine the impedance of a 
sample under test and is rewritten here only for discussion [Ref. 4, pp. 329- 
334]. 


The normalized surface impedance of the sample 1s written, 


Zour | 
ASUEE = 1/6 + ov = coth [km ~ Jaq) #/e! (2.78) 


where ¥ = |Zoyrel/eoc + ¢ being the phase angle of |Zar¢| and Zeyr¢ 


is the impedance at the surface of the sample which is to be determined. At 
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the resonance peak, the peak frequency Wr is approximately equal to the 
normal angular frequency of longitudinal vibration, uw, , Or 

Wy = Ys (2.79) 
At the half power frequencies (-3 dB down from the peak) the half-power half 


width (HPHW) can be written as 





apha =k, , (2.80) 


where kj) = k, + kya, and kn, is the spatial damping constant without the 
porous material present (i.e., it takes into account wall and driver losses). 
Therefore km, the spatial damping factor due to the sample under test can be 


written as 
kK, = ki- kk. (2.81) 
Assuming the ratio k),)/uw << 1, the surface impedance of the sampole can be 
determined once ayr gr Wo, &, and c are measured. 
3. Yaniv's Method 
Equations have been developed in transmission line theory for 
determination of the input impedance of an electrical transmission line in 
terms of the line characteristics: propagation constant, characteristic 
impedance and length of line [Ref. 5]. When the line is terminated in an 


open or short circuit, the input impedance of the line may be expressed in 


terms of the line parameters with the following equations: 


os = 2Zchar coth [yL} @ie>) 


N 


and 


me 


Zeeman a tanh || LL) (27283) 
where Zo, is the input impedance to a line terminated in an open circuit 


(Figure 2.4), Zo. is the input impedance to a_ line terminated in a short 
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circuit (Figure 2.5), Zohar 18 the characteristic impedance of the line, 
Y = jk 1s the complex propagation constant of the line andL is the length 
of the line. 


From Eq. (2.82) and Eq. (2.83) we can solve for the material's 


Characteristic impedance: 


ond ~~! wt 


2char = [Zo Zog! ae (2.84) 
and for the complex propagation constant 
tanhlyE) =. (222/2e eo (2.85) 


The high-impedance termination for an acoustical line may be realized by 
termination with a rigid wall and a low-impedance termination for an 
acoustical line may be realized by terminating the line with a column of fluid 


equal to a quarter wavelength, itself terminated in a rigid wall [Ref. 5]. 


SOovuree TER Minanhoms 





Figure 2.5 Geometry for Yaniv's Short Circuit Determination 
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By measuring the standing wave inside the tube alternatively 
terminated with an open circuit and then a short circuit, both Zo.) allceeZ as 
may be determined. With these two values, Eq. (2.84) and Eq. (2.85) may be 


solved for the desired complex propagation constant and characteristic 


impedance. 


See PROBE TUBE METHOD 

This method was first presented by Scott (1946) [Ref. 2] to allow the 
determination of the complec wave number ffrom the characteristics of an 
acoustic disturbance propagating through a fluid-saturated porous’ solid by 
probing the sound field existing in the material and extracting the wave 
number and absorption coefficient. 

To determine the complex wave number in the porous medium, which describes 
the spatial variation of the acoustic disturbance for a uniform, progressive 
plane wave, Scott writes the velocity potential in of the average motion in the 


pores as, 
ej wt (2.86) 


where k =k - ja is the complex wave number and x represents distance in the 


direction of propagation of the wave. Eq. (2.86) may be rearranged such 


that, 
~ —ax _j(wt — kx) 2207) 
? = do e7 

where a is the absorpotion coefficient and k = 27n/iX is the wave number. 


Experimentally, one measures the acoustic pressure amplitude as a_ function of 
distance in the material, along with the wave length. 
To determine the characteristic impedance of the porous material, Scott 


writes, 
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~ Befe(w)w 


Henne — i (2.83) 





where poe¢¢(w) 1s the complex frequency dependant effective density of the 
fluid in the pores. Determination of Pose: assumes of course, that the 
effective fluid density is known, for example, from Biot theory. Eq. (2.88) 
was used only to check that the measured values of k were consistent with the 
values of re from the other methods. Only the values of k measured with 


Scott's method are renorted here. 


D. MAGEREAL “PROPERPIESs iE IHdOD 

This method requires determination of five properties of the material tG | 
calculate the complex pronagation constant and characteristic imnedance of all 
material. Beranek [Ref. 4, Ch. 19] discusses this method in detail and only 
that portion ot the theory applicable to the comparison with the new technique 
will be presented. 

Beranek has shown that both the propagation constant and characteriustic 
impedance are derivable with the aid of charts fron the basic physical 
properties of the material (Ref. 4, p. 839]. 

These properties are: 

(a) the specific flow resistance 

(b) the porosity 

(c) the structure factor 

(d) the volume coefficient of elasticity of the air in the interstices 

(e) the volime coefficient of elasticity of the skeleton of the material. 
First, one must determine whether the material under investigation can be 
considered a soft blanket or a_ rigid/dense blanket. This determination is 


made by considering the ratio of the volume coefficient of elasticity of the 
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Semin tne Interstices £6 that Gf the skeleton. If this ratio is less than 


S the material is considered to be a rigid/dense blanket [Ref. 4, p. 840]. 


The fiberglass material chosen for the rigid/dense blanket criterion for 
Beranek provides an equation for the propagation constant 5 = jk, of such a 


material, namely 


~~, /ookst¥ Ry 
b = jk = Jw = Ne poKet ’ (2.89) 


where 
b = complex propagation constant (same as Yaniv's Y), 
k = complex wave number, 


kst = Structure factor (same as Blot's tortuosity, a), 


K = volume coefficient of elasticity of air in the interstices (i.e., 
the bulk modulus, and can be isothermal, adiabatic or complex), 


Ry = dynamic specific resistance per unit thickness (approximately the DC 
flow resistance), 


P9 = denisity of air, and 

Y = porosity (same as P in Biot's theory). 
By the use of Figure 2.6 and Figure 2.7 the magnitude and phase of b can be 
determined. The characteristic impedance may then be determined from 


2 =, ee 
char ay * 


2 


(2290) 


This method suffers from the fact that not all of the material properties can 
be experimentally determined and thus reasonable assumptions for their values 


must be made [Ref. 4, p. 840]. 
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Figure 2.6 Chart for Determining the Magnitude of the Propagation Constant b 
of Rigid Tiles [Figure 19.3 of Ref. 4, p. 843] 
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Figure 2.7 Chart for Determining the Phase Angle of the Propagation Constant b 
of Rigid Tiles [Figure 19.4 of Ref. 4, p. 844] 
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IIIT. EXPERIMENT AND RESULTS 


The objective of this study was to develop a new technique to measure both 
the complex wave number and characteristic impedance of a fluid-saturated 
norous solid. As with any new technique, comparison with classical methods is 
a good basis on which to judge one's theory. This section will first present 
the experimental results of the five classical techniques discussed previously, 
and then the new technique. A comparison of measured data for all methods will 


be presented in Chapter IV. 


Ae PROPERTIES OF THE POROUS SOLID 

The porous solid used in all experiments was a semi-rigid type of 
fiberglass, much like Owens-Corning Tupe 705 fiberglass, used in the insulation 
of refrigeration type equipment. This material has sufficient rigidity to 
allow it to be cut from sheets into a variety of shapes and to allow it to be 
inserted into tubes without excessive compression. This property was important 
in the flow resistance and impedance tube measurements where a snug fit was 
required. 

Discs with nominal dimensions 3-29/32 inches in diameter and 1-29/32 
inches thick were cut from three sheets of this material, weighed, and measured 
to determine the material's volume. Fram two sets of measurements, the average 
material density was found to be 43.0 + 5.0 kg/m. For comparison, 
Owens-Corning Type 705 fiberglass is reported to have a denisty of 96 kg/m? 
[Ref. 5], while rock-wool, also known commercially as Stillite, is reported to 
have a density of 80 kg/m? [Ref. 2]. 

The porosity of the fiberglass material was measured by immersing a l-inch 
cube in a graduated beaker of water and measureing the volume of displaced 
water. The fiberglass material itself was assumed to be non-absorbing. The 


resulting porosity was found to lie between 0.99 and 1.0. 
4] 


1. Flow Resistance Measurements 
Next, measurement was made of the specific DC flow resistance per unit 
length of sample. For most homogeneous, isotropic materials this’ flow 
resistance may be measured with sufficient accuracy by applying a known steady 
air flow across a_ sample and measuring the resulting pressure differential. 
The actual method used is similiar to that presented by Beranek [Ref. 4, p. 
SSui 


The quantities measured were: [Ref. 4, p. 847] 


| 


Ap = pressure differential between the two faces of the sample 


V velocity of the linear flow of alr through the sample due to Ap 
The pressure differential was measured with a standard water manometer (+ 0.1 
cm accuracy) and the velocity of linear flow was determined buy measuring the 
time required to collect a sufficient volume of water (which drew air through 
the sample by a siphon) in a graduated cylinder. The experimental apparatus 
used is shown in Figure 3.1. All tube connections were sealed with RIV Silicon 
Rubber to prevent leaks. 
The velocity of linear flow produced by Ap ls, 


VS ’ (3015 


where 


V = volume of water collected in time t, and 
A = cross-sectional area of sample. 


Then the DC flow resistance Rpc is given by 


where p is the density of water, g is the acceleration of gravity, 4h is 
the difference in water height in the two arms of the manometer, and T is the 
thickness of the sample. The units of Rpc are in Rayl/unit thickness, which 


can by written in COGS or MKS units gm/cm2*sec/cm or kg/m**sec/m, respectively. 
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In the first flow resistance experiment seven discs were placed into — 


the large sample holder, each disc having a diameter of 9.85 am and a thickness 
Oi 4.70 = cm. The samples were made as close-fitting as possible but, 
unfortunately, the die used to cut these discs had a diamter approximately 0.1] 
cm greater than that of the sample holder. While providing a tight fit along 
the tune wall, this led to a compression in the material's thickness by 
approximately 10% as they were pushed into the holder. At that time it was 
thought that the effect of this compression would not be significant enough to 
require re-machining the cuting die. Later measurements of flow resistance 
uSing a different holder, however, gave values which were a factor of 2 smaller 
than that for flow normal to the fiberglass surface. 

Measurements were conducted over a limited range of air flow velocity 
by adjusting the 3/4 inch Hoke valve on the water side of the apparatus to give 
the desired pressure differential. When steady state conditions were reached, 
the volume flow was then plotted as a function of volume flow rate. This is 
shown in Figure 3.2. Within experimental error, the variation in Ap with V was 
found to be linear up to the highest flow rate measured, which waS many orders 
of magnitude greater than expected for a sound wave. Since the sound pressures 
appearing at the surface of a sample are seldom greater than 1000 dynes/em* in 
practice (100 Pa or 0.5 inches water) it can be seen that even at higher 
pressure differenitals the linear relationship betwen Ap and V is maintained. 
The DC flow resistance had a value of 59,100 + 4,200 MKS Rayl/m using the 
uncompressed thickness. Table 2 lists the experimental data and 
Table 3 lists the (DC flow resistance results using the compressed and 
uncompressed thickness. 

It was discovered later (while attempting to reconcile these 
measurements with the DC flow resistance required to fit the resonance data) 


that the DC flow resistance of our fiberglass samples is quite anisotropic (the 
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DC flow resistance required to fit the [100] modes resonances, which has its 


particle velocity parallel to the air-fiberglass surface, was between 10,000 to 


20,000 MKS Rayl/m). 


PRESSURE DIFFERENTIAL 
(Pascal ) 


Press. Difference vs Volume Flow 
(VOLUME FLOW NORMAL TO SURFACE) 





128.38 188.8 248.8 338 .8 


VOLUME FLOW RATE 
(cc’sec) 


Figure 3.2 Ap vs Vin the Large Sample Holder 
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TABLE 


LARGE SAMPLE HOLDER DC FLOW RESISTANCE MEASUREMENT DATA 


(Flow Normal to Surface) 


Uncompressed thickness = 33.34 + 0.l cm 
Compressed thickness = 29.80 + 0.1 cm 


Sample holder diameter = 9.84 + 0.1 am 
Room temperature = 23.3 deg. Centrigrade 





Ah Vor, t V Roce 
(ec) (cc) (sec) (cc/sec) (MKS Rayl/m) 
6.40 20 e) 2038 + 4 Fico tee Ore 278.0 0.5 57 5x10 2 ee 
6.4 2 it) 2 Te6 214 5/26 

bes 2074 7 a4 278.8 Gres 

eG 2016 8.4 240 .6 58.2 

DED 2056 Cray 2360.6 58.1 

5 2008 8.9 2254 58 .8 

NS. 2046 1oe2 198 .6 S66 

a Zee TOR 2012 55:9 

As8 2108 Ieee 184.9 60.8 

5e5 2128 are 15 35 ST oe 

325 2063 14.3 144.7 60.4 

Be 2074 13.2 53.46 56.9 

Zoo 2077 Seen 114.6 54.5 

225 2062 ig<c 109.5 See 

2a 2103 19.0 1109 56.35 

ee 2098 2730 77.9 54S 

lier 2081 26.8 77.6 54.8 

ae 20¢2 26ian. Ie 5456 

eal 2086 46.9 44,4 61.8 

Ieee 2120 46 .0 Aer 65.1 

eZ 2100 44.3 47.4 OSes 

O56 yan 98.5 21.4 70.0 

O6 2116 94.4 2204 66 .8 

OxG 2161 93°.5 2362 64.7 


*Compressed Thickness Data 
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TABLE 3 


LARGE SAMPLE HOLDER DC FLOW RESISTANCE MEASUREMENT RESULTS 
(flow normal to surface) 


Compressed Thickness: 


Mean Roc = 59.1 x 103 MKS Rayl/m 


Std. Deviation = 4.2 x 103 MKS Rayl/m 
Uncompressed Thickness: 


Mean Rpc = 53.2 x 103 MKS Rayl/n 


Std. Deviation = 4.2 x 103 MKS Rayl/m 
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A smaller samnle holder was subsequently constructed and samples were cut which 
allowed the ~easurement of flow resistance with particle velocity normal to the’ 
fiberglass surface and with particle velocuty parallel to the fiberglass 
Surface. Figure 3.3 shows the results. 

Not only is the IX flow resistance anisotropic by a factor of two tom 


three but the effect of the 10% compression of the material in the large sample 


holder ts “evident. The resulting flow resistance parallel to the surface GE 


the material was found to be 15,200 + 500 MKS Rayl/m while that normal to the 
Surtace was 34,000 + 900 MKS Rayl/™. Table 4 and Table 5 give the measured | 
data. | 
2.  Viscous=Inertial Transtt lon si regucme 

The ratio of the viscous penetration depth to the characteristic pore 
size is one of the most fundamental dimensionless parameters that determines 
the nature of elastic wave propogation in a fluid-saturated porous solid. ig 
Syisc “< pore size, then the motion of the fluid in the pores is essential 
the same as if the fluid were inviscid. However, if Syjigg >> pore size then the 
fluid in the pores is essentially frozen to the solid, and sound cannot 
propagate within the fluid. Microscopic photographs of the fiberglass material 
provided an estimate of the average distance between fibers. This mean 
distance was determined to be 220 + 80 im. The crossover from where fluid 
inertia dominates to where flow resistance dominates occurs when Syjsc #* 1/2 
(mean pore dlameter). This value was estimated to be 110 +40 um. Since, the 
thermal penetration depth in air is approximately the same as the viscous 
penetration denth, the crossover from isothermal to adiabatic compression in 
the pores was expected to be nearly coincident with the viscous-inertial 
transition. Figure 3.4 shows the viscous penetration depth over the 


frequencies of interest. The transition frequency for this case was 
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BMOrOXimeate ly WN) ez. Pewee l=ckyOu OL Sttaiqht pipes of radius 4, the Xx 


Plow tosistance 1S Jiven by, 
See 
Cae a 3) 


Consideriny the '% flow resitance measured with Flow parallel co the 
ait-tiberglass surface, the equivalent radius would be 98 un. This figure 
comnares very well with the microscopic observations. A radius of 98 um results 


mia transition frequency of 500 Hz. 


peel’ SSIC "ETHODS OF TRANS IISSION LINE THEORY 
ieeeelieciatice lube Experiment and Results — 

The impedance tube technique is the most widely used to determine a 
material's characteristic imnedance, being relatively simple in both theory 
and measurement. This exneriment was the first conducted of the five classical 
methods. 

A Bruel and Kjaer Type 4002 Standing Wave apnaratus was used for these 
measurements. The tube was constructed of clear nlexiqlass which allowed 
visual inspection of the sample. The tube is 1 meter in lenath and 9.8 ~™ in 
diameter. This sets a lower limit on the usable frequency range to be 
approximately 85 Hz for the length of the tube to be greater than 1/4 
wavelength, and a higher limit of approximately 2006 Hz for operation below the 
lowest crossmode, thus ensuring only plane waves are propagated. the HP4192A 
Low Frequency Impedance Analyzer was used for signal generation and detection. 
It has a built in frequency synthesizer with a resolution of 1 mHz anda 
nhase-sensitive detector with a 0.1 Hz bandwidth to preclude noise and 
harmonics interfering with measurements. Figure 3.5 shows the experimental 
setup. Since we are primarily interested in the lowest longitudnal mode in the 
rectanqular cavity (particle velocity parallel to the surface of the 


fiberqlass), which has a frequency in iS range of 600 to 700 Hz, the impedance 


Press. Difference vs Volume Flow 
(FLOW NORMAL/PARALLEL TO SURFACE) 
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Figure 3.4 Characteristic Length Scale of a Pore Compared to the Viscous 
Penetration Depth 
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TABLE 4 
SMALL SAMPLE HOLDER DC FLOW RESISTANCE MEASUREMENT DATA AND RESULTS 
(Flow Normal to Surface) 


Thickness = 15.60 + 0.01 cm 
Cross-Sectional Area = 20.5 + 0.1 cm? 





An Vol. f V RDC 

EC ) EEC) (sec) (cc/sec) (MKS Rayl/n) 
feo t 0.) Dies, 24 Feo ee | 278-66 £ 0.5 35.1xl02 + 900 
we 2086 720 205-2 3571 

6 .0 2052 9.0 229.3 3340) 

5.9 2096 9.4 22243 342 

4.9 2082 1 ea 185.6 32-0 

50 2108 10.9 192.9 33.4 

Ae 0 2084 135 154553 Sod 

4.0 2116 1374 74 3227 

Bou 2150 18.3 7 36 3269 

50 2098 18.0 116.6 332 

2.0 2118 2en4 74.6 34 .6 

ZO Zio 29.1 74.0 34.8 

120 2116 by ac 36.6 B52 

<0 212) 55.4 38.3 B35 


103 MKS Rayl/m 


Mean kpc as S40 ex 
0.9 x 103 MKS Rayl/m 


Std Deviation = 
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TABLE 5 
SMALL SAMPLE HOLDER DC FLOW RESISTANCE MEASUREMENT DATA AND RESULTS 
(Flow Parallel to Surface) 


Taileknesse=) )/.00) 20.01 cn 
Cross-Sectional Area = 20.5 + 0.1 cm2 





Vol. t Vv Roc 
(Cc) (sec) (cc/sec) (MKS Rayl/m) 
0.1 2068 +4 Smee, | PIG.8 32 0.5 14.7x103 + 500 

2100 7.5 281.1 14.4 

2066 9.0 229.8 ae 

2120 9.1 239.2 sre 

2109 11.9 eigen 15.9 

2124 Melee 191.5 14.7 

2126 14.9 142.8 1520 

DiS 15.2 139.0 lise 

2075 202 94.0 15.6 

2108 22.9 92 <0 15.9 

2150 ADD eo 1525 

2104 39.6 53.3 14.9 


Mean Rp = eS 


2 x 102 MKS Rayl/m 
Std. Deviation = 0.5 x 


103 MKS Rayl/m 
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tube measurements were carried out at frequencies between 100-1000 Hz with 
sample thicknesses of 4, 8, and 12 inches. The different sample thicknesses ~ 
were used in order to determine how thick the material had to be for’) 
retlections from the rigid termination to be negligible, thus allowing the 
material's characteristic impedance to be determined. 

ileasurements of the reflection coefficient from 100 to 1000 Hz” forewaiae 
three samples indicated that the material could be considered infinite for 
thicknesses greater than 8 to 12 inches. The reflection coefficient, R, is 
qiven by 

R= |B/Al , (3.4) 
where the ratio B/A is given by Eq. (2.68). The reflection coefficient changes 
noticeably in the frequency region betwen 400 and 500 Hz. # £=This is thought to 
be due to the transition from isothermal to adiabatic compression of the gas in 
the pores. 

The standing wave ratio and the distance to the first pressure node 
were measured for the 12" sample at frequencies between 100 and 1000 Hz. The 
normalized characteristic impedance 7 ee Viose was calculated using the 
"Smith" chart nomogram with the SWR and phase at the first node as input 
parameters. Table 6 gives the results. 

It was observed that at frequencies below 200 Hz fluctuations occur in 
both maqnitude and phase of the impedance. The position and the amplitude of 
the pressure antinode and node were very difficult to measure accurately below 
200 Hz, and so these data are not as reliable as the data taken above 200 Hz. 
From 200 to 400 Hz, the impedance magnitude decreases at a larger rate than at 
frequencies above 500 Hz, apparent evidence of the isothermal to adiabatic 
transition occurring in the pores. This transition is seen in both maqnmameue. 


and ohase as well as the reflection cocifiieieae- The transition coincides 
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Figure 3.5 Experimental Setup for the Impedance Tube Method 
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TABLEMG 


IMPEDANCE TUBE MEASUREMENT RESULTS ON 12" OF FIBERGLASS 


Distance 4 

Freq. £o flrs ~ Zchar 
eZ SWR node (cm) lZchar! /Poe (Degrees ) 
100 oa 59.70 1.60 -44.0 
nZ5 Doo 50.70 3.74 -48 .0 
200 4.79 36 .80 2.64 —3 2a 
250 4.40 2a) o Zool —a ie 
300 4.00 24.10 ZeeS ei |. 
400 3.46 18.00 2.60 =36e0 
500 Zeos 14.45 1.95 —22009 
600 a7 11.40 1.94 —3l 
700 2.26 9.05 1.80 —2 9 
800 2-30 8225 1.75 -3 i 
900 Zs i leceill -24.4 
1000 2 ler 6.90 jer) —2 5 
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with the velocity profile changing from laminar flow in the low frequency limit 
to plug flow in the high frequency limit (see Figure 2.1). 
2. Resonance Analysis Experiment and Results 

Of all the classic mehtods, this technique in most similar to the new 
Gecinique in that parameters are extracted from the resonance curve to 
determine a sample's characteristic impedance. As previously commented, this 
method is limited to discrete frequencies and thin samples since the resonances 
of the tube become weaker as thicker samples are introduced. The upper limit 
thickness for our fiberglass was found to be about 2 inches. A l inch thick 
sample was chosen since Yaniv's method was conducted with the same sample 
thickness and since three of the five empty tube resonances were still distinct 
enough to make measurements. 

The experimental setup is very similar to. the impedance tube 
experiment. A Bruel and Kjaer Type 4002 Standing Wave apparatus was used along 
with a HP3580A Spectrum Analyzer. The spectrum analyzer was used to determine 
the resonance frequencies and half-power widths of the lowest three peaks in 
both the empty tube and in the tube filled with the 1 inch thick fiberglass 
sanple. A HP3314A Function Generator was then set to the frequency of each of 
these neak and the distance between the fiberqlass surface and the pressure 
maximum closest to the speaker was measured. This 1s needed to calculate the 
effective tube length, as an infinite impedance loudspeaker 1s assumed in the 
theory. The surface impedance of the sample was then determined from these 
measurements using Eq. (2.78). Figure 3.6 shows the experimental setun and 
Table 7 gives the resulting values of normalized impedance values. 

The lowest two resonances (center frequencies of 253 Hz and 400 Hz, 
respectively) gave plausible values of surface impedance, however the highest 


mode with the absorbing material present (center frequency of 538 - 544 Hz) had 
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Figure 3.6 Equipment Setup For Resonance Analysis Method 
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TABLE 7 


RESONANCE ANALYSIS SURFACE IMPEDANCE RESULTS 
(Impedance Tube and Rectangular Cavity) 
(1 Inch Thick Fiberglass Sample) 


Surface Impedance Results in Impedance Tube 


Center 

Frequency 

(damped ) HPEW ~ Zourf 

(Hz) (Hz) IZ sure |/Po¢ (deqrees) 
253 23 Zo leaves. 
400 18 Zeek —88 .9 
543 41 enh +66 .4 


Surface Impedance Results in Rectangular Cavity 


Center 

Frequency 

(damped) HPFW s Zourf 

(Hz) (Hz) IZ sure |/ Poe (degrees) 
841 126 lst -65.4 


This point was discarded, as discussed in text. 
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0.0078 
0.0013 


=O 00> 


Kpp/4n 
0.071 


a higher Q than in the empty tube and was therefore discarded. The erroneous” 
data for the third resonance peak is attributed to the significant broadening 
of its resonance curve and the resultant overlapping with adjacent modes. 

This resonance analysis method was applied to the rectangular cavity 
with and without al inch thick strip of fiberglass while exciting the [010] 
mode to verify the results obtained in the tube. The [010] mode was chosen | 
because for this mode the particle velocity is normal to the fiberglass 
surface, aS in the impedance tube. For these meaSurements in the rectangular 
cavity the source was placed at grid point (0,0) (Figure 3.11) and the receiver 
at grid. point (20,0). The resonance curve was first measured in the empty 
Cavity; it had a center frequency of 873 Hz. The normalized magnitude of the 
surface impedance obtained agreed closely with the value obtained by the 
impedance tube for frequencies between 800 and 900 Hz. However a Significant 
phase discrepency was evident. The reason for the discrepency 1S not known. 
Perhaps 1t 1S a result of uSing a polnt source located ina corner. Table 7 
lists the results. 

As with other impedance tube measurements, this method relies upon a 
plane propagating wave, where incidence on the porous material is assumed 
normal to the surface of the sample under test. With a point source located in 
one corner of the cavity used, departure from plane wave propagation can be 
expected with increasing sample thickness as the quality factor of the 
resonance degrades. The extracted impedances are surface values. In een 
geometries in which this method was used the magnitude of the surface impedance 
decreased with increasing frequency (which was expected). The phase in all 
measurements indicates that the impedance remains largely reactive, and that 
the gas in the pores acts like a_ stiffness. Use of samples of greater 


thickness would only suppress the resonance curves more. 


ao 


3. Yaniv's “lethod, Experiment and Results 

The third and last technique which uses transmission line theory is the 
method pronosed and experimentally verified by Yaniv [Ref. 5]. The experiment 
was conducted using the impedance tuhe of the previous two experiments, with 
several modifications. Yaniv's method requires measuring the surface impedance 
with the sample backed by a rigid termination and by an acoustic short-circuit, 
obtained by moving the rigid termination such that a 1/4 A air gap exists 
between the rear surface of the material and the rigid termination. 

The rigid end cap used in the two previous experiments was renlaced 
with a 2 inch thick disc of plexiglass with a 3/8 inch diameter hollow aluminum 
tube threaded into the side away from the driver. A small amount of high 
vacuum grease was used to seal the edges of the disc at the tube walls while 
fim) allowing for positioning. The front face of the disc had a recess cut 
out of it into which a FET microphone mounted at the end of a 3/16 inch hollow 
aluminum rod could fit snugly and flush so as to provide a smooth surface when 
the disc was placed against the back end of the sample. The 3/16 inch diameter 
rod fit inside the 3/8 inch diameter rod, providing support and allowing the 
microphone to be positioned anywhere between the backside of the porous solid 
and the rigid termination. Electrical connection to a power supply switchbox 
was also made by wires which passed through the small hollow rod. = This 
configuration allowed the disc and microphone to serve as the rigid termination 
for the open Circuit impedance measurement and also allowed accurate 
positioning of the termination to provide the required one-quarter wavelength 
air gap between the porous solid and the rigid termination at each frequency. 
Figure 3.7 shows the movable termination and Figure 3.8 shows the equipment set 
setup. 

A 1 inch thick disc of fiberglass was inserted into the imnedance tube 


Such that the distance betwen the driver and the sample was at least a quarter 
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Figure 3.7 The Rigid Sliding Termination (Yaniv'S Method) 


wavelength. The lowest usable frequency was therefore approximately 250 Hz, 
the upper limit remaining approximately 2000 Hz to preclude the lowest 
crossmode. Measurements were made of the first minimum and maximum pressure 
amplitudes in front of the sample and of the distance between the sample 
Surface and the first pressure minimum in fron of the sample. For each 
frequency, these measurements were performed under conditions of ooen and short 
Circuit terminations. A short circuit termination at the back side of the 
fiberglass disc was obtained by adjusting the position of the plexiglass plug 
until a null output was observed from the microphone, which was held flush to 
the fiberglass disc. It was found that, within experimental error, the result- 
ing distance between the porous disc and the plexiglass plug was within 1% 
of the 1/4 wavelength determination made in the empty tube at the frequencies 
used in the experiment. The standing wave ratio and the reflection coefficient 


phase angle were calculated using Eq. (2.67) and Eq. (2.70), respectively. The 
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Figure 3.8 Equipment Setup for Yaniv's Method 
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open and short circuit impedances were found using the "Smith" chart. The — 
Characteristic impedance was then found by use of Eq. (2.84). 
Determination of the complex propagation constant, y = jk, required 


finding the roots of Eq. (2.85) such that, 


tanh[yL] - [Zeo/Zog) 1/2 = 0 . (3.5) 
The roots were found using a numerical technique known as Ward's method. The 
details are given in Appendix A [Ref. 19]. The results of this experiment are 
listed in Table 8. 

As for the previous methods, a decrease in impedance magnitude with 
increasing frequency was observed, with the transition from isothermal to 
adiabatic compression between 400 to 700 Hz in this 1 inch thick sample. The 
complex wave number also shows this transition. The data and results for 400 


HZ appear to anomalous. The reason is not known. 


C. PROBE TUBE METHOD, EXPERIMENT AND RESULTS 

The probe tube measuremnets were made by setting up a plane wave in the 
Bruel and Kjaer Type 4002 Standing Wave apparatus filled partially with 
fiberglass and exploring the sound dield by means of a probe tube extension to 
the impedance tube microphone. The slope of the decay of the pressure 
amplitude with distance provided a measurement of the attenuation constant 4a. 
The change in phase of the oscillatory wave with distance provided the wave 
number k. The complex wave number is then found from k =k - ja. 

The equipment setup is shown in Figure 3.9. Three specimens, of 
thickness 4, 8, and 12 inches, were chosen for measurement and comparison with 
the impedance tube results. Fach sample thickness was made up of the required 
number of fiberglass discs, 9.84 an diameter and 4.76 cm thick, inserted into 


the tube opposite the driver, taking care to avoid compression of each disc 
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while trying to minimize any air gaps between disc surfaces. Each disc hadi 
0.71 cm diameter hole cut in the center for a close fit around the probe tube 
without any deformation of the fiberglass. 

The probe tube was manufactured out of steel tubing with the same outside 
and inside diameter of the existing open ended probe tube on the standing wave 
apparatus. A copper cap was welded onto the end of the extension next toa 
3/16 inch diameter hole to allow the measurement of the sound field. The 
extension attached to the existing probe tube with a plastic sleeve. This is 
shown in Figure 3.10. 

The experiment was conducted at frequencies of 125, 250, 500, and 1000 Hz 
for each of the three sample thicknesses, but only the 12 inch thick sample 
was long enough so that the wave reflected from the rigid end of the tube was 
sufficiently attenuated to be negligible compared to the incident wave. The 
pressure amplitude at the rigid termination was measured to be 20 dB less than 
that at the sample boundary. Determination of the attenuation constant was 
made by plotting pressure versus distance on semi-log paper and extracting the 
Slope. The phase of the pressure was measured with a SR-350 Dual Channel 
Lock~In Amplifier. The wave number was determined from the change of phase 
with distance into the fiberglass. Table 9 displays the experimental results 
in terms of k and ak. 

The impedance magnitude decreases with frequency while the phase indicates 
that the material becomes more resistive with increasing frequency. Only the 
absorption coefficient shows any sign of the transition between isothermal and 


adiabatic compressions in the pore volume. 


D. MATERIAL PROPERTY METHOD, EXPERIMENT AND RESULTS 
All the methods discussed so far involve measuring certain characteristics 


of a progressive plane wave in the porous material or its proximity. Beranek 


65 


KLKUSUL HP S4C0A 
DSS 5020 RMS 


DEGe AL VATMETE ro 
STORAd & 


OSLTLLD SCOPE 










SR S30 
LOCK - TN 
AMPLIFLER 










WP SSI4A 
FUNCTION 
GENERATOR 






ITWACO 1201 
Low NOLSE 


PREAMPLIFIER 

















HP ALFA 
POWER AMPLIFIER 


> 


HP SSOOB 


FREQUENCY — 
COUNTER, 






HP SSBO0A 
SPECTRUM 
ANALY ZER 








OK, 


One® 


PROBE TUBE 
MTECROPHONE 


Figure 3.9 Equipment Setup for Probe Tube Experiment 


66 


has shown that the complex wave number and characteristic impedance can also be 
derived from the basic physical properties of the material with the aid of 
charts. These properties are: 

a) specific flow resistance per unit length, R)(* Rp) 

b) porosity, Y 

c) structure factor, kg, (tortuosity) 


ad) volume coefficient of elasticity of air in the interstices K (Bulk 
modulus) 


e) volume coefficient of elasticity of the skeleton of the material. 
Measurements of the specific flow resistance and porosity have already been 
discussed. The remaining three material properties are difficult if not 
impossible to measure, and will, based on a reasonable Knowledge of the 
materlal, be approximated. 

The structure factor kgy, known as the tortuosity in Biot theory, is 
essentially a path length correction factor. Experimentally, it is found to 
depend most strongly on the porosity of the porous’ solid [Ref. 11]. In the 
limit porosity = 1, the tortuosity also equals 1 for typical porous merials, 
and certainly for fiberglass. As the porosity of the samples used in this 
research was found to be very close to l, a value of 1 for the tortuosity was 
assumed in all calculations. 

The volume coefficient of elasticity of air in the interstices, K, 1s an 
important term; if isothermal conditions can be assumed then K = 1 x 10® dynes/ 
cm?; if adiabatic conditions exist then K =1.4 x 10° dynes/cm’. For most 
fibrouS acoustic materials however, K will, in some frequency region, be 
neither isothermal nor adiabatic but will be something in between and will be 
complex valued and frequency dependent, as is our case. 

Finally, the volume coefficient of elasticity of the skeleton, Keke}, 
although not directly involved in the calculation of k and 7 eae must be 


known or approximated since it is the ratio of K/Kepe) that determined whether 
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a material is treated as a soft blanket or a rigid/dence blanket. Based on the 
appearance, texture, and high flow resistivity of the fiberglass it was assumed 
that the material should be treated as a rigid/dense blanket. (Ref. 4, pp. 
839-840] 

Once the properties of the material were characterized, Eq. (2.89) could be 
used to directly calculate or Figure 2.6 and Figure 2.7 used to graphically 
determine the magnitude and phase of the complex propagation constant b = 3k as 
a function of frequency, from which the complex wave number is determined. 
Once these were determined the characteristic impedance could be calculated by 
use of Eq. (2.90). The values determined from the graphs_ for k and the 
calculated values of epee assuming isothermal conditions in the pore volume 
are listed in Table 10. The same quantities, assuming adiabatic conditions, 
are listed in Table ll. 

As expected, the impedance magnitude decreases with increasing frequency. 
The impedance phase calculated for both isothermal and adiabatic conditions is 
the same, a consequence of the bulk modulus being a_ real quantity in the two 
limiting cases. If isothermal conditions were to exist, then the impedance 
phase would approach an angle of -45 degrees (pure diffusion of sound 
pressure). Since the calculated phase values remain nearer to the diffusive 
limit, the assumption that isothermal conditions exist in the pore volume would 
be most nearly correct. 

Over the frequency range considered, the isothermal values for k and a are 
greater than the adiabatic values, which 1S consistent with the fact that the 


1sothermal sound speed is less than the adiabatic sound speed. 


E. THE NEW TECHNIQUE, EXPERIMENT AND RESULTS 
The plane wave resonator chosen to conduct all experiments was constructed 


out of plexiglass with interior dimensions as follows: 
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Length (x-direction) = 25.4-2 0.1 aca 

fidth (y-direction). — 196 0ies 

Height (z-direction) s——9 2392-0. bene 
This particular cavity was constructed for a laboratory experiment on standing 
waves in a rectangular room. The minimum wall thickness was 1.27 a™m (1/2"). 
The cavity was mounted on three leqs with the source hole located in the floor 
at one corner to allow the excitation of all modes. The top cover of the Gavia 
is removable to allow insertion of fiberglass samples. When installed, the ton 
cover seals the cavity and slides to allow the receiver position to be varied. 
A small micronhone receptacle mounted on aé=é slidebar installed in the 
removable cover allows the pressure to be measured throughout the cavity with a 
probe tube. All sliding surfaces were constructed with very close tolerances 
and coated with high vacuum grease to reduce leakage and to provide a tight 
fit when the cavity cover was installed. 

The receiver slide bar assembly consisted of three portions: the slide bar, 
the microphone receptacle, and the receiver microphone. The slide bar fit into 
a groove machined into the interior surface of the top cover and fit flush with 
this surface. A threaded hole was cut into the slidebar to allow mounting the 
microphone receptacle. The microphone receptacle consisted of a hollow 
aluminum cylinder with a thumb screw to set and hold the height of a brass 
probe tube assembly. The probe tube assembly allowed insertion of a General 
Radio precision microphone flush with the smaller bore probe tube extension 
which could be inserted into the cavity interior without disrupting the sound 
field. 

The bottom of the cavity was marked with grid lines, 21 in the length 
direction, 15 in the width direction, spaced 1/2" apart, producing 336 


individual Grid Pperats. The in-phase and quadrature components of pressure 
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were measured at each grid point for the mode and sample selected. Figure 3.11 
shows the location of the source and the grid numbering scheme that was applied 
Berthie interior Eloor of the cavity. 

Ten fiberqlass samples were cut from a fiberglass sheet, with thicknesses 
Beer) Erom 0.4e con (3/16 wnch) to 4.76 cm (1 7/8 inch) the maximum thickness 
heingqg equal to the full thickness of the sheet. The samples were cut in the 
form of rectangular strips so as to fit Snugly along the side wall opnosite the 
source. Figure 3.11 shows the location. The strips extended from qrid point 
mmi>) tO grid point (20715). The samoles were cut such that the surface 
normal to the flow velocity in the original flow resistance experiment would 
also be normal to the particle velocity expected when the [010] mode was 
excited (i.e., the "width" direction in Figure 3.11) and parallel to the 
Marticle velocity when the [100] mode was excited (the "length" direction). 
Table 12 lists the samples used and their thicknesses. 

Three experiments were conducted: 

1. The mode structure in the cavity was measured over the frequency range 
of interest for each sample. The effect of introducing the fiberglass 
samples on the strength and shape of the resonances was noted. Of 
particular interest were the [100] and [010] modes, which were expected 
to be relatively distinct. 

2. The DC flow resistance of the air-filled fiberglass was then extracted 
from the above measurements using the dispersion relation developed 
in Chapter 2 for the allowed modes (Eq. (2.51)), following which the 
complex wave number and characteristic impedance were calculated. 

3. The pressure field of the [100} and [010] modes was measured for 
varying sample thicknesses and compared to that predicted by theory. 

The equipment setup used to conduct all three of these experiments is _ shown 


ipeeagure 3.12. 
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Figure 3.11 Rectangular Cavity Grid Map 


1. Modal Frequencies Up to 2000 Hz for Various Fiberglass Thicknesses 


The first of the three experiments was conducted to determine the modes 
present in the empty cavity and the effects of increasing sample thickness on 
the shape and strength of each resonant mode. From mode theory, seven modes 
were expected to exist in the empty cavity for frequencies less than 2000 Hz. 
Of these seven modes, the two lowest were the [100] and [010] modes, both of 
which occur below 1000 Hz. All seven modes were observed, with the two lowest 
modes being strong and relatively distinct compared to the others. This is 
shown in Figure 3.13. Table 13 lists the theoretically calculated and 
experimentally measured peak frequencies of the modes present in the cavity up 


£© 2000 4Hzs 


Te 





There was some concern that as thicker samples were introduced the 
resonance peaks would overlap, making the measurement of the half-power widths 
unreliable. This condition set the upper limit on usable sample thickness to 1 
7/8" for the [100] mode (sample #10) and 11/16" for the [010] mode (sample 
po). 

For all observed modes, the center frequency down-shifted as the sample 
thickness was increased, while the Q decreased. The [100] mode has an empty 
cavity peak frequency and O of 680.2 Hz and 147.9, respectively. Inserting the 
1 7/8" (sample #10) fiberglass strip down-shifted its center frequency and 
lowered the 0 to 570.0 Hz and 3.2, respectively, for a total frequency change 
of 110.2 Hz. The [010] mode has a center frequency of 870.2 Hz and a OQ of 57.6 
in the empty cavity, but when the 11/16" (sample #5) fiberglass strip was 
inserted the peak frequency only down-shifted 17 Hz to 853.2 Hz with a O of 
Z\al. For both mode measurements the receiver probe tube remained flush with 
the ceiling of the cavity at grid point (20,0). See Figures 3.14 and 3.15 

Several observations from this experiment were made: 

(1) The Q and the center frequency of both modes decreased with increasing 
sample thickness. 

(2) Compared to the [010] mode, the [100] mode, with particle velocity 
parallel to the fiberglass surface, allowed thicker samples’ to be 
introduced before the resonance curve spread such that the 0 could no 
longer be reliably determined. 

(3) For all sample thicknesses, the Q of the resonance curves’ was 
larger with the receiver positioned above the source (grid point (0,0)) 
than in the opposite corner (grid point (20,0)). This was assumed 
to be due to near field effects of the source on the receiver. 
Consequently, only the resonance data taken with the receiver at grid 


point (20,0) is reported here. 
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Figure 3.12 Experimental Setup for the New Method 
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(4) Overlap of the (110) mode with the [010] mode was evident for 
fiberglass thickness of 1/2" (sample #4) and greater. Overlap of 
the [100] mode with the [010] mode was not’ significant until a 
thickness of 1 1/4" (sample #8) of fiberglass was introduced. See 


Figures 3.14 and 3.15. 


2. Measured f, and _ OQ of the [100] and [010] Modes for Various Fiberglass 


Thicknesses 

The second experiment consisted of measuring the center frequency and 0 
of the [100] and [010] modes as a function of sample thickness. The receiver 
was placed in the ceiling of the cavity at grid point (20,0) at the opposite 
end of the long wall from the source. This location was chosen to minimize 
direct acouStic cross-talk between the source and receiver. For the [100] 
mode, each of the 10 samples was inserted and the center frequency and O were 
measured from the resonance curve. This sample procedure was repeated 
for the [010] mode except only samples 1 through 5 were’ used. The center 
frequency and half-power frequencies were mesured to within a tenth of a Hertz 
by reading the frequency of the 1 to 1.5 MHz Local Oscillator output of the 
HP3580A Spectrum Analyzer on a HP5300B Frequency Counter. Setting the scale to 
1 dB/div on the Spectrum Analyzer allowed the center frequency and half~-power 
full width of the resonance curves to be measured. Tables 14 and 15 list the 
experimentally determined values of f, and Q for the [100] and [010] modes. 

Following the experimental measurements, theoretical values of f, and 
Q for the [100] and [010] modes as a function of sample thickness were 
calculated fram the complex roots of the dispersion relation (Eq. (2.51)) for 
various values of DC flow resistance using a numerical analysis technique Known 
as Ward's method [Ref. 19, pp. 169-170]. A brief description of this method is 


given in Appendix A. However, before a comparison between theory = and 
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Figure 3.14 Mode Structure of Cavity With a 11/16" Fiberglass Strip 
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experiment could be made, and the best fittinq value of DC flow resistance 
determined, 1t was necessary to correct one or the other for empty cavity 
losses. Possible sources of empty cavity losses are considered in the next 
subsection. 
3. Sources of Loss in the Empty Cavity 
In deriving the dispersion relation for the new resonance technique it 
was assumed that the only losses which occur in the plane wave resonator are 
due to viscous and thermal losses within the fluid saturated porous solid. 
Although these are the largest sources of loss inthe cavity, other los 
mechanisms must be accounted for. To account for these, several assumptions 
were made: 
(1) the losses in the empty cavity remain the same even though a 
portion of the volume and wall surface is lost due to the introduction 
of a fiberglass sample (less than 10%), 
(2) the shift of resonance frequency due to the introduction of a_ fiber- 
glass sample does not effect the Q of the empty cavity and, 
(3) each loss mechanism can be treated as acting independently. 


Generally, stated, it was assumed that, 


Total losses in the = losses in the + losses due to 
fiberglass filled empty cavity only to fiber- 
cavity glass 


which can be written in terms of quality factors as 


ote (3.6) 
where Onyp is the Q of the cavity due to all loss mechanisms, Qp is the empty 
cavity Q and Op is the Q due to viscous and thermal losses within the 
fiberglass. Eq. (3.6) was used to compare the measured values of Qppp to the 
calculated values of Op for the [100] and [010] modes by setting the values of 
Or equal to those experimentally determined for the empty cavity. 
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Sample # 


Empty 


10 


[100] MODE CENTER FREQUENCY AND Q_ VALUES 


Thickness 
(cm) 


0.0 


0.48 


ee, 
lag 
Zece 
2.94 


3.18 


Cc 
Hz) 





680 .2 


669.5 


667.5 


658.6 


652.1] 


646.3 


623.6 


614.1 


52° 


D/126 


50.0 


26 


TABLE 14 
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67S ae 


663.9 


oG 3 


648.8 


640.1 


628.9 


587.0 


567.4 


532 


506 .6 


495.4 


+ Q 


ol 





674.1 


668 .1 


664 SZ 


652.2 


05055 


647.6 


O27 a 


624 .3 


672.6 


5262 


34.1 


27 ame 


2072 


9.8 


lou 


6.1 


4.9 
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Sample # 


Empty 


(010) MODE CENTER FREQUENCY AND Q_ VALUES 


Thickness 
(em) 


eu 
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i> 





TABLE 15 


Hz) 





863 .0 


854.5 


852.0 


846.4 


841.1 


834.4 
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eee) 


ot 





DG 


43.5 


Orr 


Boe2 


26.8 


Zi 


+ 0.9 


| 


An attempt to identify the major source of loss in the empty cavity was 
made, as this in turn determined the value of Op to be used in correcting for 
losses other than in the fiberglass. The loss mechanisms in the empty cavity 
were assumed to be due to 3 mechanisms: 

-- transmission of energy through the plexiglass walls of the cavity, 
-- viscous and thermal losses at the interior surfaces of the cavity, | 


-- leakage of energy due to incomplete sealing of the top cover and sliding 
section of the cavity and slide bar assembly. 


The transmission of energy through the walls was of concern mainly since _ the 
cavity was not designed specifically for this experiment. For normaly} 
incidence, the transmission loss through the walls of the cavity was calculated | 
to be 32 dB for the minimum wall thickness (1.27 am), [Ref. 17, p. 128]. Based 
on the maximum wall thickness (2.54 om) the lowest expected coincidence 
frequency for a free plexiglass plate was calculated to be 1282 Hz. For a 
plate fixed about its edges this frequency is higher. Therefore, transmission 


by coupling to flexural waves in the walls can be ignored. The following 


equation was used to estimate the QO resulting from each source of loss: 


ee 2™ Estored (3.7) 
~  Elost/cycle 

where 
Ee eee (3.8) 


and «€ is. the average energy density, and V is' the interior volume of the 
cavity. The energy lost per cycle can be written 


. AC it 
Elost/cycle 7 ell = ’ (3.9) 


where A is the surface area, c is the speed of sound in the fluid, Ty the 
power transmission coefficient and f is the frequency. Inserting Eq. (3us 


and (3.9) into Eq. (3.7) results in 
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8 nV 


Otrans = \TTA : (3210)) 


For loss due to transmission throught the cavity walls, assuming 1000 
Hz as the worst case (A = 0.343 m), the value of OQrrang was calculated to be 
2700. It was concluded that transmission through the walls was not 
significant. 

Next considered were the effects of the viscous and thermal losses 
occurring at the cavity surfaces. Again Eq. (3.7) was used to determine Ooyrf, 
Piet beto-eq Cemaining the same as inEq. (3.8). The energy lost per 
cycle due to viscous effects can be written approximately as 

Elost/cycle = ™ %visc (3.11) 
where dyjicc 18 the viscous penetration depth in air. The energy lost per 
cycle due to thermal effect can be written approximately as 

Elost/eycle = SA (y¥-1) Sth (oie) 

where y is the ratio of specific heats Cp/Cy (= 1.4) Forme aie anc o-;, 1s 
the thermal penetration depths in alr. Since the viscous and thermal 
penetration depths in air are approximately equal, an expression for the 
energy lost per cycle due to both viscous and thermal effects is 

Elost/cycle = “ %Visc Y > (3.13) 
In other words, 

Syisc * (11) Sth = Y8visc : (3.14) 
Sueostituting Eq. (3.8) and (3.13) into Eq. (3.7) results in 


0 21V 
surf ~ ASviscY * 


R 


(3215) 


Choosing a frequency of 100 Hz for the worst case (greatest viscous and thermal 
penetration depths), with the air temperature in the cavity assumed to be 20°C, 
the resulting Ogyr¢ was 450. Therefore, viscous and thermal wall losses are 
negligible. 
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Lastly, the loss due to leakage past cavity slideing surfaces was 
considered. Since this loss could not be directly calculated, a measure of 
its significance was estimated by pressing down on the top cover of the cavity 
at the wall edges while retaking the [100] mode resonance curve data for each 
sample. The empty cavity O increased from 147 to 252, with a 0.9 Hz upshift in 
the center frequency. From this experiment, a value of Q for leakage was 
estimated to be approximately 280. This indicated that leakage was the major 
source Of loss in the empty cavity by at least a factor of two. The leakage 
problem could not be corrected without making major changes in the cavity 
design. No attempt was made to account for the upshift in the center frequency 
that occurred while pressing down on the cavity top (the frequency upshift was 
on the order of 1 Hz for each sample, while the frequency downshift due to 
introducing the fiberglass was on the average of 10 to 12 Hz drop per sample). 

4. Determination of Rpc from f. and QO 

For the [100] mode, the iterative process used to extract the [C flow 
resistance of the fiberglass from the dispersion relation (Eq. (2.51)) was 
started by choosing Rpc equal to 59,200 MKS Rayl/m (the mean flow resistance 
determined in the large sample holder with particle velocity normal to the 
air-fiberglass surface). After entering the fluid and cavity parameters (see 
Appendix A for parameters and values used), the roots of the dispersion 
relation for i were then determined for each sample using the measured values 
of f,. and Q as an initial starting point. The resulting roots for each sample 
were numerically computed until the sum of the magnitudes of the real and 
imaginery parts of the dispersion relation (Eq. 2.51) was less than l x 107 &% 
By this means, calculated values of Op vs fe were obtained for a particular 
value of Roc. 

The calculated and measured values of Qvs f, were compared in two 
ways. Fiqure 3.16 shows Oppr vs £,- Here the calculated values of Q, have 
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been corrected for empty cavity losses using Eq. (3.6). Figure 3.17 shows Op 
vs f.. Here the measured values of Opp have been corrected for empty cavity 


losses using 


ea 1 


OF or 

If the calculated curve plotted above the measured curve, this 
indicated too high a value of Roc and a lower value of Ry was chosen. If the 
curve plotted below the measured curve, this indicated too low a value of Ro 
and a larger value of Rpc was chosen. This process was continued until the 
measured curve was bounded between two values of Roc» Once bounded, smaller 
increments were chosen for Rnc and the process was continued until a best fit 
(graphically) was obtained between the calculated and measured curves 
of Op Or Op vs f,- For the [100] mode, the best fit value for Ry was 
15,000 MKS Rayl/m, which was four times less than that measured experimentally. 

Two possible causes were investigated to explain this discrepency. The 
first possible cause is the 10% compression of the samples observed in the 
Original DC flow resistance experiment (large sample holder). Assuming the 
total pressure drop across the test sample was independent of the compressed 
thickness, the extracted value of DC flow resistance per unit length should be 
corrected from 59,200 MKS Rayl/m to 53,000 MKS Rayl/m. The second possible 
cause 1s that the flow resistance of the fiberglas is not isotropic, and that 
1t 1S greater for flow perpendicular to the sheet’ faces. This might be 
expected, as there 1S visual evidence that the fiberglass sheets have been 
compressed. In the first flow resistance measurement, the fiberglass sample 
was oriented so that the fiberglass sheet surface was normal to. the flow 
velocity. However, for the [100] mode, the particle velocity 1s parallel to 
the fiberglass sheet surface. Flow resistance tests were again conducted, 
uSing a much smaller sample holder, with samples cut so that the anisotropy in 
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Rpc could be measured. As indicated earlier in this chapter, the flow 
resistance normal to the sheet surface now measured approximately three times 
greater than that parallel to the sheet surface. The value measured with flow 
parallel to the surface was 15,200 MKS Rayl/m. With this value of Rpg, the © 
roots of the dispersion relation were recalculated, Opgr and Op calculated and 
plotted. These values fell within 8% of measured values for Oppr vS f. and | 
within 12% for Op vs £, for all samples. 

The method of extracting Rnc from measured values of f. and Q was then 
applied to the [010] mode, which has particle velocity normal to the sheet 
surface. Starting with a value of 34,000 MKS Rayl/m for Roc (the value 
measured in the small sample holder), values of QO vs f, were calculated for 
various values of Roc and were compared with measured values of O vs fu. The 
measured values of O vs f, are bounded by the calculated curves for 
Rpc = 5,000 MKS Rayl/m and Rye = 20,000 MKS Rayl/m, and are consistent with a 
value of 15,000 MKS Rayl/m for Roc, the same as for the [100] mode. This is 
shown in Figures 3.18 and 3.19. The reason for the discrepency between these 
bounding values and the measured value of Roc for flow normal to the fiberglass 
Sheet face 1s not known. It may be due to mode overlap of both the [110] and 
[100] modes on to the [010] mode. Figures 3.20 and 3.21 show plots of the 
resonance center frequency for the [100] and [010] modes, respectively, asa 
function of sample thickness for the same values of Rp aS in Figures 3.18 
and oer 

Having determined that the best fit value of the DC flow resistance in 
the porous material was approximately 15,000 MKS Rayl/m, the complex wave 
number and characteristic impedance were calcualted using Eq. (2.56) and Eae 
(2.57) and the roots of the dispersion relation, (Eq. (2.51)). Table 16 lists 
the results for the [100] mode, while Table 17 lists those for the [010] mode. 
The data for both modes is graphically presented in Figures 3.22 through 3.26. 
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5. Sensitivity of Roots of the Dispersion Relation to Variation in 
Parameters 


In order to determine the sensitivity of the roots of the dispersion 
relation to variations in fluid and porous material physical properties, the 
roots were again calculated for the [100] mode with each of the input 
parameters changed by a small amount, ome at a time, and the characteristic 
impedance and complex wave number were recalculated. 

It was found that there was less than a 1% change in k and ee 
when porosity was decreased to 0.95 (5% change), when bulk fluid sound’ speed 
was changed to 343 m/s (0.6% change or 2 m/s), and when bulk fluid density was 
changed to 1.23 Kg/m? (1.6% change or 0.02 Kg/m3). However, an 11% change in k 
and Fee ay was found when the DC flow resistance was varied by + 5,000 MKS 
Rayl/m about 15,000 MKS Rayl/m. Because the uncertainty in the measured data 
is so much greater, the sensitivity of the best fit value of Rp, and hence k 
and are to slight uncertainties in other physical properties is negligible. 

6. Modal Pressure Field in a Rectangular Cavity with Porous Acoustic 

Material 

To theoretically determine the [100] and _ [010] modal pressure 
distribution in the cavity for each fiberglass sample requires the values of 
ce and Ko, be determined from the roots of the dispersion relation for "i 
for the best fit value of Ro; Then, uSing Eq. (2.63) for the bulk fluid and 
Eq. (2.64) for the porous material, the normalized magnitude and phase of the 
pressure at each grid point can be determined. See Appendix 8B for a more 
detailed discussion. For the [100] mode, Figure 3.27 shows graphically the 
results for sample #1 (3/16") and Figure 3.28 shows the results for sample #10 
igln7 7/3"). For the [010] mode, Figure 3.31 shows graphically the results for 


sample #1 (3/16") and Figure 3.32 shows the results for sample #5 (11/16"). 
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Figure 3.22 \Zchar|/P9¢ vs f£ For New Method ([100] and [010] Modes) 
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Figure 3.25 a vs f For New Method ([{100] and [010] Modes) 
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TABLE 16 
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Zchar AND k RESULTS FOR [100] MODE 


(Roc = 15,200 MKS Rayl/m) 
Ec k a _ Ue ee Soh 
Sample # (Hz) (Rad/m) (Np/m) lZchar|/PoC (degrees) (m/s) 
I, 669.4 +0.1 2].54 a .o5 1.48 ete eZ 
2 665.4 21.47 14.47 1.48 See 194.8 
3 O5/./ Ao 4235 1.49 = all 850 
4 649.8 Deal? 14.21 io =o erk bozo 
> oor. 0 Zao e275 Wee oa) Nearer 
6 O25 .3 20.74 ore be53 OOo 189.4 
y Or) al Z0e 0! LS o2 1.54 -36.6 eek 
8 S0le. 5 ZU 208 Sen oo = Ore 185.4 
e 587.4 20.19 12.66 heoc =) 6 Juste. 
10 Boo. 20.03 ine o6 1.61 -34.4 N Fish als: 
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TABLE. “17 


Zchar AND kK RESULTS FOR [010] MODE 


(Roc = 15,000 MKS Rayl/m) 
fc k }. ee ere 
Sample # (Hz) (Rad/m) (Np/m) |Zchar|/Po¢ (degrees) 
it 865.0 =) O26 P07 LO -34.8 
Z 862.9 ZZ 15.63 1.30 -3428 
3 Sou 6. 25.64 eo) 1.31 =34..3 
+ 854.3 ot 15.46 lr -34.8 
5 848 .3 25.47 15.29 1.32 =34 7) 
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These plots present the mode's modal pressure field in phasor notation. nae 
maqnitude is represented by the size of the arrow and the phase angle by the 
angular rotation of the arrow in deqrees from the horizontal (i.e., zero 
deqrees of phase angle is a phasor pointing horizontally to the right). 

Experimentally, the pressure field was measured by inserting each 
sample into the cavity and exciting the desired mode at its center frequency 
for that sample. Ry programming a HP-300 series computer to take data from the 
SR 530 Dual Channel Lock-In Amplifier (in the X, Y mode), the in-nhase and 
quadrature components of pressure were measured at each grid point. The 
normalized maqnitude and phase of the pressure was then plotted. For the [190] 
mode, Figure 3.29 displays the results for sample #1 (3/16") and Figure 3.30 
displays the results for sample #10 (1 7/8"). The results for the [010] mode 
are displayed in Figures 3.33 and 3.34 for sample #1 (3/16") and sample #5 
(11/16"), respectively. 

In comparing the theoretical and experimental pressure field plots 
several observations can be made: 

(1) The effects of cavity leakage at the slidinq surfaces is evident 
when the experimental plot (Figure 3.29) is compared to the theoretical 
plot (Figure 3.27) of the [100] mode for sample #1 (3/16") by noting 
that the pressure phasors along the boundaries not in contact with the 
fiberqlass strip in Figure 3.29 exibit a small increase in phase angle 
(i1.e., they are not horizontal as they would be if no losses existed at 
the boundaries) which is not present in Figure 3.27. This effect may 
also be seen for the [010] mode in Figure 3.31 and Figure 3.33. 

(2) The effects of the porous material on the pressure distribution in 
the bulk fluid can be seen to "grow" with increasing sample thickness by 


observing the relative increase in the phase angle away from the 


102 





14 atdures 40J vOoTzNqtAySIq eANssaead Sepow [NOT] TeoTzEeAOeUTL Lz°E eaNbtTy 


(----— 








— ee Oe Oe ei 


ee 


( WOWIXHW OL A3IZIIBWYON) 
( 1# etdues ) ,9T/E ‘apoy [00 


Y 
aINSSetlg Tepow Teotzeiceauy, 


a 





» 5 3 é a 6 
30.ynNes Woes NOILDSYIG HLGIM ---- 


103 


OL# eTdwes AzoJ UOT INGTAISTG eaNnsseag epow [NNT] TeoTIeACeUT 8z7°€ sanbty 


ee S3IYNOS WOYS NOTLOSYIC HLONSI 








| 
| 
| | 
| 
| 


—s of Cd t-— — ee, n,n, ne a 
__ oo ) -  Or  OO —> > (= — ee 
-_—_———oo oS OO oO) Oo —~ a" @ — om -—- a n,n, ee, ee 
Arepunog sseTtsueqty-dutv 
-_—— + —-—)—— — ob ¢ / oe —— —— —— ee (——__* 


| 
| 
: 
| 
| 
| 


L L L L L L 4 4 t 1 d d 
. 7 [ameter tot ot t tT t <t f 
(amie ete | «6FlCOUTlCUECE 


| 
| 
| 
| 
| 
{ 
eee eee. ts ly 


———— A a SS 


Bl 2 9 S b € 2 t Q e 
( WOWIXEW OL D3ZITWWYON) ° 
( OT# eTdues ) yB/d-T ‘SPOW [00T] “C1 .08! 
PTeTyI sunssstg Tepow Teo Isesios9uy, et 





104 


See MOL fo BOG. Gt thet 
35ynNos WOus NOILO3SYIC HLGOIM ----> 





I# aTdures 20J YOTANGTAASIG 


---- 39 


Arepunog ssel[dslaqtyq-d 


ynos WONIS NOILIINIA 


8ANSSeAd epow [ONT] Tequowtrzsdxq 67°F 


HLONIY 


oanbt 4 





= 
at nd gy = o eee 
a ~ Ol OO 
lg lg ly ~ - ee lee 
a ee aed, 
a a a ek: Gay Ee 
ee ln Ly sy ~ o SHO eae i eee eee eee eee eee 
ee a e SO ee 
eee ly Ly Ly ~ - ee 
ee ly lg le _ e (es ee eee 
ne lg Ly Ly oe ees ee 
el) Ug Wy Ly Ce eee ee 
lg Wg Ly, Ly e ae ea 
td gas OS Oo 
——— 2 ee lg Ll Se ee a ee ee 
ee Seen yes es ay =. pone ae es oe ee 
2 
YON 
( T# tdi wes -) ‘O1/C ‘3 3p0 a 
PLety ainssadg Tepow pe 


z ¢ » & S§ 2 6% & Gt tt 
3DyNOS WONs NOILDIYIMC HIGIM ----> 


105 


Ol# aTdues AOJ UOTINYIAASTG sANsseig apow [OOT] Tequewirzedxy oE'E oanbry 


(———— Senos HOYs NOILIS8IC HLONAI 


sane 


Azyepunog sSseTssaqty-dUuyy 
a Ee? 4 es i gue = 


‘ole ky 
- 








8 OS a 
er NN ONDA Rw 
\ 
, 
| 


c— =e et oe 


20unoc 


i—_- b 8 t£ 9° § F&F ¢ 2 
( WIHIXEM OL D3ZITENUON) 
( o1# etdwes ) «B/2-1 ‘apow (00 T] 
PLeTYy FdnssstIg Tepoy Paeinseoy 


06ST bi OSI 


t 8 b Ql 


S08nosS WONd NOILO3SeTad HLGIM 


I 27 ¢« Fv §¢ 4 


> 


Oo 
2° 
x = 
el ot 
3 


106 


———————_- °° °° °° 


I# eTdwes 10} YOTINGTAISIG aansseig epow [OTO] TeoTJeso09uL, TEE eanbt 4 


(<<-- 3JDaNOS WONS NOTLISYICG HLON|I 


Sey 








Axrepunog sseTsZseqty-uty 








—_—_ --_ -_~ -——_—_->— =) 2 Ly SD» 'D ty 
rt re ee 


—) —_— —> —> —- — —> — — — — — —S —> — ——b — ad se — — 


ee eee ee eee 
meee —a—e—eeeia———Ca ais eg 
—_ eC eer a a a OCR Cet 
— oon nr tO a eR 


rr rm re nner rn ne nner a pF 











| 6 8 

( WNWIXYUW OL GSZITIBWYON) 
( t# otdwes ) ,91/€ ‘poy [Oto] 
PT29TY FINnssodlg Tepow TROT asOdu], 





€ *© S SS ¢ Me TSeeGi atest 
JDYNOS NOMS NOPLISstG eh aa 


107 


G# aTdwes AO} UOT|ANGTAASTG eaANssarag spow [OTTO] TeotqeAoeuL zZEe°E aanbty 


(---= 304yNO0S HOY NOILDSYIG HLONZNI 





ed ay 





Axsepunog sselssueqty-upty 


I eee |. 


| — —— o— — — +— —. (= — — — —— — — — —— -_ (——— — <a eee 
| -—_- --—_ -—- + 8 —  ——  —  — CU 
TT TTT rr rr rr a 
rr fri ff 
) 
eens ee a BFR ine c 
Coeee ci cle Ct NCO c! | st | ef). zt 81 ) ay c 8 d 9 5 » € 2 Q opt 
( WAWIXYW OL GAZITWWYON) p 
( S# etTdwes ) ,9I/Tt ‘apow [oto] ey) ve 
PLPTY eINSSeatIg TeRpow TeOT IssO9U], s 


€ » Ss 3 é 
alana nOaseNOT basta 4iediM -—--—-> 


108 


Td atdues aoyz uOTANYTAASTG eansselag epow [OIO] Tequewtazedxq EE°¢g oanbrg 


——__ == JDYNOS WONd NOILIBYIG HLON]AI 


Axepunog sseTSuaqyTy-duty 





———) <——__y <<  , ,, E—EEeeee 
— — > > > > SP ns» _, —— —___, ——., ——_, 
— —-_—o—o—-——iaseeo— avr" a seo—_—n se ss ees snr ee se Sen ees es eee Ose 
eT a _e eee eee ees se em ee ss ees eee ees eee asi 


—) —- oa —) —- —) a=) ) =} <= —) =) oS ob oad sy —) oa) ~y oma) ah 


—— ¢ < a ee 





") 


9ounos 


i @ 

( WIXWW O1 AIZITIYWYON 
( 1# atdues ) 191/E. "a poLl [O10 
PLT°TYy euNnssadg Tepoy pounseoay 


+ SS” & 2) @ 398 eect ot 
3DYxNOS WOYS NOTLISYTA HLGIM ---~-> 


g 
© | 
ss % 
y > 
+ *o 
2 


Co 


air-fiberqlass boundary from the theoretical plot of the 3/16" (sample 

#1) fiberglass strip compared to the theoretical plot of the 1 7/8" 

(sample #10) fiberglass strip. This can be seen in Figures 3.27 and 

3.28, respectively. This effect is not seen in the theoretical plots of 

the [010] mode for sample #1 (Figure 3.31) and for sample #5 (Figure 

S52) 

(3) The effects of mode overlap become more evident as sample thickness is 
increased by comparing the measured and calculated modal pressure fields 
for the same sample thickness. For example, compare Figure 3.29 (sample 
#1) to Figure 3.30 (sample #10) for the [100] mode and Figure 3.33 
(sample #1) to Figure 3.34 (sample #5) for the [010] mode. 

(4) The pressure distribution at the air-fibergqlass boundary shows’ the 
effect of the relative size and sign of kj, and k5, by noting that for 
the thickest samples the pressure phasor's angle increases as- the 
interface 1S approached. For the [100] mode commare the theoretical nlot 
for sample #1 (Figure 3.27) to the theoretical plot for sample #10. 
(Figure 3.28) and for the [010] mode compare the theoretical plot for 
sample #1 (Figure 3.31) to the theoretical plot for sample #5 (Fiqure 
B2 32). 

The primary reason for comparing the calculated and measured pressure 
distributions was to verify which mode was excited and to check its purity. The 
determination of whether the expected mode had been excited was made Dy 
comparing the calculated and measured pressure field plots for the thinnest 
fiberglass strip (3/16"). For the [100] mode the nodal plane runs normal to the 
Peyeeeecdirection from grid point (10,0) to grid point (10,15). Travelling 
across the nodal plane, a 180 degree phase shift in the pressure phasor is 


observed, with the pressure maqnitude reaching a maximum at the cavity 
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boundaries normal to the length direction. Likewise, for the [010] mode the 
nodal nlane is normal to the width direction falling half-way between the 
boundaries parallel to the length direction. The pressure phase also changes 
180 deqrees in crossing this nodal plane, with its magnitude reaching a maximum 
at the fiberglass boundary and opnosite wall. 

The purity of these modes is retained only with the thinner samples and 
becomes contaminated as thicker samples are introduced. This is clearly seen by 


eomparing Figures 3.29 and 3.30. 


Jy 


IV. COMPARISON OF RESULTS AND CONCLUSION 


A. COMPARISON OF THE COMPLEX WAVE NUMBER 





The complex wave number may be written in terms of its real and imaginery 
parts as the wave number k and the absorption coefficient a , respectively. 
Figure 4.1 shows the wave number k as a function of frequency for each method. | 
Values of k determined by the new method fall between those found by Yaniv's” 
method (the upper bound) and by Beranek's Material Properties method in the 
isothermal limit (the lower bound). For example, at 600 Hz the value of the wave 
number determined by the new method was within 6% of the value determined by 
Yaniv's method. The values of k from the probe tube experiment are about twice 
those determined by the new method. This is to be expected, since the particle 
velocity in the probe tube experiment was normal to the air-fiberglass boundary, 
i.e., waS in the direction of higher measured flow’ resistance. (Note that 
although the [010] mode also has particle velocity normal to the fiberglass 
surface, the best fit value of Rpc (15,000 MKS Rayl/m) was used to calculate both 
the complex wave number and characteristic impedance). 

The abosorption coefficient as a function of frequency is displayed in Figure 
4.2 for each method. The values determined by the new method fall about those 
obtained by Beranek's Material Properties method (in the isothermal limit). For 
example, at 600 Hz the value of the absorption coefficient determined from the 
new method falls with 17% of the value determined by Yaniv's method. Again, the 
values determined using the probe tube method were about twice those measured by 
the new method. The cause of the dip in the values of the absorption coefficient 
determined by Yaniv's method 1s not known. 

Figure 4.3 displays the phase speed in the porous material as a_ function of 
frequency for each method. Values determined by the new method are bounded by 
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those found by Yaniv's method and by Beranek's Material Properties method (in 

the isothermal limit). For example, at 600 Hz the phase speed determined by the 
new method 1s within 7% of the value determined by Yaniv's method. The values | 
of phase speed determined by the probe tube method are approximately one-half of 


those found by the new method for the same reason discussed above. 


Be. COMPARISON OF THE CHARACTERISTIC IMPEDANCE 

The normalized characteristic impedance as a function of frequency for all 
methods is displayed in Figure 4.4. Values determined by the new method are 
bounded by those from the adiabatic and the isothermal material properties 
method. Values from the new method plot almost directly on top of the Biot | 
theory in the diffusion limit. For example, at 600 Hz the value determined by 
the new method fell within 26% of that determined by Yaniv's method and within 
20% of the impedance tube value. 

Figure 4.5 displays the impedance phase as a function of frequency for each 
method. The values determined by the new method fall about the values 
determined by the material properties method. The low frequency values (< 200 
Hz) determined by the impedance tube method are thought to be in error due to 
the difficulties in measuring the distance to the first node in that method. It 
1S not known why the values of the phase extracted using the new method from the 
[100] mode data increase for decreasing frequency, since in the low frequency 


limit, the phase should approach -45 degrees (pure diffusion of sound pressure). 


C... CONGESTION 

The ability to determine a material's DC flow resistance from parameters of 
the resonance curve in a plane wave resonator was verified. Additionally, the 
values of both the complex wave number and characteristic impedance determined 
from the new method compared well with classical methods; the values obtained by 


the various methods agree with each ego to within 26%. 
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APPENDIX A 


METHOD OF SOLUTION OF THE DISPERSION RELATION AND ITS IMPLEMENTATION 


Hoe witteitecdmised tO Lind the roots CF Ehe dispersion relation (Ea. (2.51)) 
was Ward's method [Ref. 19, pp. 169-170]. This metnod is not the only method 
available to find the roots of a complex transcendental equation but is easy to 


understand and to implement, as no deriviatives are needed. 


1. WARD'S METHOD 
If we let the polynomial f(z) be a function of the complex variable 


meee + iy such that 


Bean =_ulx;y) > 1Vixy) + eee lee) 
femean be Shewn that fitz) has a zero a = b+ id if and only if u(b,d) = v(b,d) = 
Dig This pair of simultaneous equations in u(x,y) and v(x,y) may be 


conveniently solved in two ways. First, one can seek to minimize 
u*(x,y) + v2(x,y) [Aal.2) 


or, secondly, on can minimize 





w(x,y) = |u(x,y)| + |v(x,y)]. (ee .3) 
Since no derivitives are employed, this method seems well suited for finding 
multiple roots. 

Ward's technique is based on the knowledge that in every neighborhood of a 
point (x,y) such that w(x,y) #0, there exists a point w(x',y') such that w(x', 
y') < w(x,y). At a given step in the process, which is lterative, the current 
value of w(x,y) 1S compared with the vlues w(xth,y), w(x-h,y), w(x,yth) and 
w(x,y-h) until a smaller value is found (1f possible). The coordinates yielding 


this smaller value of x and y are used in the next iterative step. Should these 


four points yield no decrease in w= |u| + |v|, the increment h_ is decreased 
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by some amount and the process is repeated. As with other methods which can sho 
; 
convergence within certain limits, there is no quarantee of convergence 


with this method. 


2. IMPLEMENTING WARD'S METHOD 

In order to apply Ward's method to find the roots of the dispersion 
relation, it must be separated into its real and imaginary parts. To do so note 
that the parameters K1> and Oy, and as i) are functions of w , they may be 
separated into their respective real and imaginary components (these three terms 
also have complex parameters buried within them which are frequency dependent). 
Although tedious, it is necessary with this method. 

Once each parameter has been reduced to its simplest complex form (real 
and imaginary portions), these can be inserted back into the dispersion equation 
and by use of complex arithmetic, the dispersion relation may be recast in the 
form 

Re {Dispersion Relation} + jim{Disperson Relation} = 0 
and Ward's Method of Eq. (A.1.3) may be now implemented in a computer algorithm. 

The numerical analysis program written using Ward's Method required the 
following inputs for each sample thickness: 

(1) Measured Quantities 

a) ££. (measured ingthesgecayi ty) 
b) OQ (measured in the cavity) 
(2) Cavity Parameters 
a) Length of the cavity, L. (in x-direction) 
b) Width of porous material, dj. (sample thickness) 
c) Width of remaining bulk fluid, dj- 
d) Mode number in the long direction. {n=0 for the [010] mode, 
n = 1 for the [100] mode} 


Al 


(3) Fluid and Porous Solid Parameters 
a) Adiabatic fluid sound speed, cg. 


b) Fluid bulk density, p¢- 


(@ 23°C, 345.2 ms) 


aaezo-G, 1.2) Ko/m?) 


Seite ar teu cabio Of Specific heat of fluid, y. (4 20°C, 1.402) 


qd) Prandtl number, Pr. 
e) Porosity, P. (= 1.0) 
Program Parameters 

a) Angular frequency increment, Aw. 
b) Damping frequency increment, A8. 
c) An ending criterion for w(w,8). 


Initial Estimate of Roc 
a) If 


(i) 


no a priori 


knowledge of material's 


Came OM Err 7 33) 


flow resistance, then 


assume very hiogh or low value as the starting point. 


b) Perform a CC flow 


for Rm: 


resistance experiment to extract a mean value 


Once these quantities are entered, the program calculates its first w(w,8) based 


on the initial inputs of f, and Q. 


four incremental values of wlw + dw, 8 
initial w(w,8) with 
w( wt Aw, B) 
w( uw- Aw, B) 
w(w, B+ AB) 
w(w, B-A8) 


Which ever is the lowest, those coordinates in 


values and the process is 


Criterion. The output converts the roots of 


22 


continued until w(w,8) 


This value of w(w,8) will be compared to the 


+ Ag), that is, it compares the 


w,® Space are chosen as the new 
is less than the ending 
and 0O. 


w and 86 to f, 


This process is continued for each sample thickness, and the resulting 
values can be plotted either as Qvs f£, or as f£,vs sample thickness and 
compared with the respective curve after loss corrections to Q are made. An 
evaluation is made as to whether Roc 1s too high or too low, the value of Roc 
changed and the process repeated for this new value of Ry for all samples, 
This lterative process 1S completed when the roots are found for a value of Roe 
which provides the best fit curve through the measured data. With this value of 
Rpc and the roots for each sample thickness determined, the characteristic 
impedance and complex wave number may be calculated by use of Eq. (2.56) and Eq. 


(Ze57) 
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APPENDIX B 


DETERMINATION OF THE THEORETICAL MODAL PRESSURE DISTRIBUTION IN A PLANE 


WAVE RESONATOR WITH AN ACOUSTICALLY ABSORBING MATERIAL 


The equation that governs the pressure distributuion in the bulk fluid 
mzeqion 1) is: 


ae ~~ ~ soot 
$5) = —jupecoslk,x] ¢10 cos[k,,(Zz-d, )] el , Sol ol) 


The equation that governs the pressure distribution in the porous material 


(region 2) is, after substituting Eq. (2.62) for $99: 


es _ cos[k, dy} - iat 
= eee ee icosen | cos [k,, (z+d, )] e (ies 
222 


With the value of Rpc known, along with the roots of the dispersion relation for 
a given sample thickness, the pressure distribution in region 1 and region 2 
can be calculated. 

In the experimental measurement of the modal pressure field, the quadrature 
component of the pressure was zeroed at the reference position of (0,0,1) to set 
a reference point such that both theory and experiment could be compared. The 
0,0 indicates the grid point (0,0) and the 1 indicates that the receiver probe 
tube was located in the ceiling of the cavity. To set the reference at x = 0, z 
=eeymecuen that the phase of the pressure distribution at this point is Zero, 
requires that 

$19 = 1 /-AIN(-w/8) . (B.1.3) 

To calculate the theoretical modal pressure distribution, start at the 
reference location of x = 0 and z= dy and incrementally move away from this 
position in increments of Ax and Az such that the theoretical grid ponies match 
those existing in the actual cavity. 
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